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Abstract: We have previously introduced a hybrid strong/weak coupling model for jet quench¬ 
ing in heavy ion collisions in which we describe the production and fragmentation of jets at weak 
coupling, using Pythia, and describe the rate at which each parton in the jet shower loses en¬ 
ergy as it propagates through the strongly coupled plasma, dEjdx, using an expression computed 
holographically at strong coupling. The model has a single free parameter that we fit to a single 
experimental measurement. We then confront our model with experimental data on many other 
jet observables, focusing in this paper on boson-jet observables, finding that it provides a good 
description of present jet data. Next, we provide the predictions of our hybrid model for many 
measurements to come, including those for inclusive jet, dijet, photon-jet and Z-jet observables in 
heavy ion collisions with energy ^/s = 5.02 ATeV coming soon at the LHC. As the statistical un¬ 
certainties on near-future measurements of photon-jet observables are expected to be much smaller 
than those in present data, with about an order of magnitude more photon-jet events expected, 
predictions for these observables are particularly important. We find that most of our pre- and post¬ 
dictions do not depend sensitively on the form we choose for the rate of energy loss dE/dx of the 
partons in the shower. This gives our predictions considerable robustness. To better discriminate 
between possible forms for the rate of energy loss, though, we must turn to intrajet observables. 
Here, we focus on ratios of fragmentation functions. We close with a suggestion for a particular ra¬ 
tio, between the fragmentation functions of inclusive and associated jets with the same kinematics 
in the same collisions, which is particularly sensitive to the x- and E’-dependence of dE/dx, and 
hence may be used to learn which mechanism of parton energy loss best describes the quenching 
of jets. 
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1 Introduction 

The LHC has ushered in a new era in the exploration of the properties of matter under extreme 
conditions. By colliding Pb ions at center of mass energies in the multi-TeV regime, the LHC has 
provided us with droplets of the hottest matter ever produced in the laboratory and a diverse suite 
of copious high energy probes with which to explore the microscopic properties of the strongly 
coupled, liquid, quark-gluon plasma (QGP) discovered at RHIC [1-4]. For jet probes in particular, 
the results of the successful LHC Run 1 have shown that many of the properties of these funda¬ 
mental QCD objects are substantially modified when the jets are produced in Pb-Pb collisions as 
compared to when they are produced in p-p collisions [5-20]. These modifications are results of 
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the final state interaction between the jets and the droplets of hot QCD matter formed in heavy ion 
collisions. As the principal underlying effect is the energy loss suffered by each of the components 
of the jet showers on their way out of the hot matter, the various modifications of jet properties 
observed in heavy ion collisions are referred to, in sum, as jet quenching. The phenomenon of jet 
quenching was first discovered without reconstructing individual jets via the strong reduction in 
the number of intermediate-px hadrons in heavy ion collisions at RHIC [21, 22]. Precisely because 
varied modifications of jet properties can now be well measured, this suite of probes has the poten¬ 
tial to provide us with unique and important information about the properties of QGP and about the 
interaction between energetic partons and this strongly coupled liquid. The imminent start of the 
LHC heavy ion Run 2, which will increase both the center of mass energy and luminosity for heavy 
ion collisions and hence will substantially increase the production rates of jets and all other hard 
probes including high energy photons and Z-bosons, makes many much more quantitative analyses 
of the striking phenomena observed in the first LHC run imminent. Turning this opportunity into 
precise extractions of QGP properties from the experimental data to come requires a diverse suite 
of theoretical tools. In this paper, we shall present substantial advances in the development of one 
such tool, introduced in Ref. [23]. 

It is the discovery that the QGP of QCD is a strongly coupled liquid, with intense collective 
phenomena, no apparent quasiparticle structure, and a very rich phenomenology that makes it of 
such interest. But, at the same time, it makes the theoretical description of its properties and 
dynamics much more challenging, since many of the perturbative tools available to describe weakly 
coupled hard QCD processes become inapplicable in the strongly coupled liquid plasma produced 
in the range of temperatures accessible to colliders. Fortunately, in recent years gauge/gravity 
duality has emerged as a tool with which to analyze the dynamics of the strongly coupled plasma 
in various non-Abelian gauge theories. Although it is not yet known whether QCD itself has a dual 
gravitational description, calculations done using these new holographic methods (where one maps 
a question in the gauge theory plasma onto a calculation done in its dual gravitational description) 
have yielded many qualitative insights into the properties and dynamics of the QGP in QCD, as 
produced and probed in heavy ion collisions. (For a review see Ref. [24].) 

Although we now have a toolkit that includes perturbative methods valid for weakly coupled 
hard processes and holographic calculations that are valuable guides to strongly coupled dynam¬ 
ics, the description of the interactions between jets and the plasma remains challenging because 
physics at both hard and soft scales enters in central ways. Jets in the 100 GeV range are produced 
in hard processes with large momentum transfers which must be described by perturbative QCD. 
Furthermore, the energetic, virtual, partons produced in a hard scattering immediately experience 
a rapid sequence of branching processes which reduce their large initial virtuality, ultimately down 
to the hadronization scale. The evolution of this branching process is governed by the DGLAP 
equations and is described well by perturbation theory, since the initial virtualities of the partons 
in the shower that results from the branching are large. Nevertheless, when the jet fragments in¬ 
teract with the plasma created in Pb-Pb collisions, soft exchanges of energy and momentum with 
the plasma constituents are governed by physics at scales of order the temperature of the medium, 
where a strongly coupled description is mandated. And, it is these soft interactions that are critical 
to understanding the modification of jets produced in Pb-Pb collisions caused by their interaction 
with the plasma. This multi-scale nature of the physics of jet probes and their modification means 
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that neither approaches that seek to treat the both the jets and their interaction with the plasma 
perturbatively [25-59] nor those that seek to address the entire dynamics of the quenching of ener¬ 
getic probes at strong coupling via gauge/gravity duality [60-78] can capture the physics at all the 
relevant scales correctly. 

To face this challenge, in Ref. [23] we introduced a phenomenological model specifically 
designed to treat the different dynamics arising on different scales differently, using a perturbative 
description where it is appropriate and deploying insights from holographic calculations where soft 
momentum transfers and strongly coupled physics enter. While there have been other attempts to 
combine results obtained from weak and strong coupling [79-85], our approach is distinct since 
it focusses on using two different calculational frameworks at the different energy scales where 
each is appropriate. In a nutshell, our model combines a perturbative description of the creation 
and evolution of jets as implemented in modern Monte Carlo codes, in particular in Pythia [86], 
with a prescription for the mechanism by which, and consequently the rate dEjdx at which, each 
parton in the Pythia shower loses energy that we infer from the holographic computation [77] 
of how energetic light quarks lose energy in the strongly coupled plasma of a gauge theory with a 
dual gravitational description. This hybrid model contains a single parameter which controls the 
stopping distance of fast partons in the QCD plasma. We fitted this parameter to data in Ref. [23], 
and then confronted this hybrid strong/weak coupling model with data on several single-jet and 
dijet observables, obtaining very satisfactory agreement. Nonetheless, our analysis also showed 
that this same set of data could be described almost as well by different assumptions for the energy 
loss rate dE/dx of particles in plasma. This indicates that the observables considered in Ref. [23] 
are not very sensitive to the microscopic dynamics governing the interaction between an energetic 
probe and the strongly coupled plasma. 

In this paper, we continue our exploration of the predictions of the hybrid model that we 
introduced in Ref. [23], and whose construction and main aspects we review briefly in Section 2.1. 
In Section 2.2 and Appendix A, we improve some aspects of the implementation of the model: we 
incorporate the effects of the motion of the hydrodynamic fluid on the rate of energy loss, and we 
embed the jets in our model in the viscous hydrodynamic simulations of Refs. [87, 88], rather than 
the inviscid hydrodynamic simulations from Ref. [89] that we used in Ref. [23]. We discuss the 
implications of these improvements for the extraction of the value of the single parameter in our 
model in Sections 2.3 and 2.4. 

After laying this groundwork, in Section 3 we use our model to analyze photon-jet and Z-jet 
correlations (generically, boson-jet correlations) in heavy ion collisions. These observables provide 
a different and complementary set of data with which to confront our model; doing so is the first 
of our principal goals. After specifying the details of how we generate and quench our samples 
of boson-jet events in Section 3.1, in Section 3.2 we compare the predictions of our model to 
presently available photon-jet data from LHC heavy ion collisions. We find that our model is just 
as successful here as it was at describing inclusive jet observables and dijet observables in Ref. [23] 
and Appendix B.^ These many successes give us confidence in our approach. 

In Section 3 and Appendix D, we find that, within current experimental and theoretical uncer- 

Tn Appendix B, we update our model calculations of the observables from Ref. [23] to include the improved model 
implementation introduced in Section 2.2. 
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tainties, two different control hypotheses for the form of the rate of energy loss dEjdx yield almost 
as good descriptions of the single-jet, dijet, and photon-jet observables that we have investigated as 
we obtain from our hybrid strong/weak coupling model. Similarly, the gamma-jet data we consider 
in this work can also be described well by the perturbative treatment of Ref. [47]. 

Therefore, these observables are not particularly effective at discriminating between rather 
different assumptions about the interactions between the components of the jets and the strongly 
coupled medium. (This observation was first made for the specific case of the dijet asymmetry ob¬ 
servable in Ref. [90].) However, we have observed various small differences among the predictions 
of our hybrid strong/weak coupling model and those of our two control models, which opens the 
possibility that the higher energy and higher statistics data from the coming LHC heavy ion Run 2 
will yield more precise measurements that can be used to constrain the microscopic dynamics of 
how jets interact with the plasma and, ultimately, the microscopic properties of the plasma itself. 
With these goals in mind, and to facilitate direct comparisons to data that is coming soon, in Sec¬ 
tion 3.2 we provide predictions for photon-jet observables in Pb-Pb collisions at ^/s = 5.02 ATeV 
and in Section 3.3 we do the same for Z-jet observables. We include the predictions of our model 
for a number of single-jet and dijet observables at ^/s — 5.02 ATeV in Appendix C. We hope that 
near-future high-statistics measurements will make the myriad predictions of our hybrid model for 
heavy ion collisions at ^/s = 5.02 ATeV among the most incisive results of this paper. 

With the goal of finding observables that are more discriminating, in the sense that they are 
more sensitive to different assumptions about the rate of energy loss dEjdx and hence are better 
able to give us information about the microscopic dynamics behind jet quenching, in Section 4 we 
turn to exploring more differential observables. Before adding further physics (and further param¬ 
eters) to our hybrid model, we cannot use it to describe the modification of the angular shapes of 
jets. We can, however, analyze fragmentation functions. In Section 4.1 we study the fragmentation 
functions of the jets produced in association with either a photon or a Z-boson, providing predic¬ 
tions for the ratios of these fragmentation functions in Pb-Pb collisions with ^/s = 5.02 ATeV to 
the same fragmentation functions in p-p collisions. These predictions are qualitative at best, how¬ 
ever, because they are sensitive to differences between the hadronization of jets in Pb-Pb and p-p 
collisions that are not under theoretical control. And, these observables turn out not to be partic¬ 
ularly discriminating. In Section 4.2 we use the tools that we have constructed to identify a new 
observable that is discriminating. This observable is the ratio between the fragmentation functions 
of inclusive jets in Pb-Pb collisions and the fragmentation functions of the associated (lower en- 
ergy) jets from a dijet pair. Much of the uncertainty associated with hadronization should cancel 
in this ratio, since we nowhere need p-p reference data and since we compare the fragmentation 
functions of associated jets and inclusive jets with the same range of energies, in the same Pb-Pb 
collisions. Furthermore, by explicit calculation we find a clear and substantial separation between 
the predictions of our hybrid strong/weak coupling model for this ratio and the predictions that we 
obtain from our two control models. And, we find a good qualitative understanding of why this 
separation occurs, which is to say of why this ratio is such a discriminating observable. We look 
forward to measurements of this observable in the coming LHC heavy ion Run 2 as these data 
should tell us much about the dynamics of how jets interact with the strongly coupled plasma. 

In Section 5 we discuss the main conclusions of this study and look ahead. 
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2 Description of the Model and its Implementation 
2.1 The Hybrid Model Approach 

In this Section, we briefly review the main aspects of the hybrid model that we introduced in our 
previous publication [23]. Our model is a response to the challenge of addressing the wide range 
of energy scales involved in the production of energetic quarks, gluons, photons and Z-bosons 
in hard processes followed by the interaction of the quarks and gluons with the strongly coupled 
matter formed in ultrarelativistic heavy ion collisions. Our hybrid model treats the weakly coupled 
short distance dynamics involved in the creation and hard evolution of jets perturbatively while 
at the same time using insights from holographic calculations that assume strong coupling in the 
treatment of the soft long distance interactions of the jet with the strongly coupled fluid created 
in the collision. The former processes occur at momentum scales set by the large virtuality of 
the elementary hard partonic interaction, which we shall denote Q. The latter processes involve 
momenta at or close to the typical scales that characterize the strongly coupled fluid, which we 
take to be of order the temperature T. For the energetic processes involved in jet production in Pb- 
Pb collisions at the LHC, Q ^ T. We use this separation of scales to justify our explicit separation 
of the treatment of the relevant dynamics for each regime. We are, however, certainly not doing a 
systematic expansion order by order in powers of the ratio of these scales. Instead, we construct a 
model in which the qualitatively different physics at the widely separated scales is each treated in 
the most appropriate language, in so doing introducing a single model parameter that at the present 
stage we must fit to data. 

In our model, we assume that the soft in-medium processes cannot alter either the short dis¬ 
tance production of hard partons in elementary partonic collisions or the hard branching processes 
by which the large initial virtuality relaxes. In p-p collisions in vacuum, those two processes are 
both described well by perturbative QCD; in particular, the radiative branching processes are in¬ 
corporated into evolution equations which are in turn the basis of the different high energy Monte 
Carlo event generators like Pythia [86], which we shall employ. The probabilistic implementa¬ 
tion of this evolution describes the formation of a shower of partons initiated from the initial hard 
process, with splitting probabilities dictated by perturbative QCD. In the environment created in 
a heavy ion collision, as this parton shower develops its constituents are continuously exchanging 
momentum with the matter produced in the collision. In our model, we assume that the splitting 
probabilities are not modified by these soft exchanges. This does not mean, though, that the shower 
remains unchanged. As a consequence of their soft exchanges of momenta with the medium, the 
particles in the shower lose energy and pick up momentum transverse to their original direction as 
they propagate through the strongly coupled liquid. In this publication, as in Ref. [23], we focus 
entirely on the loss of energy of the partons in the shower, as we can do so in a model with only a 
single free parameter. In a future publication we shall introduce the transverse momentum kicks, 
at the expense of adding at least one more parameter. Here, as in Ref. [23], we shall only de¬ 
scribe experimental observables that are sensitive to the loss of energy by the partons in the shower 
and insensitive to the transverse momentum kicks that they also experience. Since the momentum 
transfers between the shower partons and the medium are not large, the physics of energy loss must 
be described at strong coupling. After assigning a life-time to each parton in the shower according 
to a formation time argument [23, 91], we know the points in spacetime where the parton on each 
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branch of the shower is formed, and splits. We model the energy loss of each parton in the shower 
as a continuous process, supplementing the in-medium evolution with an explicit energy loss rate 
dE/dx that models the strongly coupled dynamics of parton energy loss. We do not track what 
becomes of the energy lost by each parton in the shower, implicitly assuming that the lost energy 
is incorporated into the strongly coupled fluid, ultimately becoming soft hadrons with momenta of 
order T that we do not model. The form that we assume for the rate of energy loss dE jdx therefore 
fully encodes all the strongly coupled in-medium dynamics incorporated in our model. 

In our model, we explore the consequences of an energy loss rate dEjdx whose form is that 
appropriate for the rate of energy loss of an energetic massless quark (excitation in the fundamental 
color representation) traversing a slab of plasma with temperature T and thickness x in the strongly 
coupled plasma of A/^ = 4 supersymmetric Yang-Mills (SYM) theory [77], 


dE 


dx 


strongly coupled 


4 1 

^in 2 /- ’ 

^ ®stop 


1 eH^ 

“ 2^r4/3 ’ 


( 2 . 1 ) 


obtained via the gauge/gravity duality. Here, E\^ is the initial energy that the massless quark has 
before it enters the plasma, E{x) is the energy that it has after traversing the slab of thickness x, 
and Xstop is the stopping distance of the high energy excitation — the smallest slab thickness that 
results in the energetic excitation losing all of its energy within the slab of plasma. In A/" = 4 
SYM theory, the dimensionless constant /^sc appearing in the expression for Xgtop is determined 
explicitly in terms of the T Hooft coupling A and is /^sc = 1.05A^/^ [67, 72, 75]. The premise of 
our hybrid model is that the form of dE j dx in the strongly coupled quark-gluon plasma of QCD 
is the same as in (2.1); we shall see that this hypothesis is uncontradicted by many and varied 
sets of data. However, there is no reason at all to expect that the relationship between and 
A should be the same in QCD and A/^ = 4 SYM theory, as the strongly coupled plasmas of the 
two theories have different, and differently many, microscopic degrees of freedom. Furthermore, 
there are ambiguities in the definition of jets in A/^ = 4 SYM theory: since hard processes in this 
theory do not produce jets [65, 92], different theoretical calculations have been developed in which 
highly energetic colored excitations are formed in different ways - no one of which is preferred 
over others as a model for jets in QCD since none is model for jet production in QCD. And, the 
proportionality constant between and Xgtop can depend on details of the particular way 

in which a highly energetic colored excitation is formed. For both these reasons, and as discussed 
in more detail in Ref. [23], in our model we will assume that any differences between dEjdx 
in the strongly coupled plasmas of QCD and J\f — SYM theory can be absorbed in the value 
of /^sc 5 which we therefore take as a free parameter whose value must be fixed by fitting to data. 
We will refer to the form (2.1) for dEjdx as strongly coupled energy loss. Our hybrid model 
constitutes applying this prescription for energy loss branch-by-branch to the partons in a shower 
that described in vacuum by Pythia. We shall specify the implementation of our hybrid model 
fully in subsequent subsections. 

In order to have some other benchmarks against which to compare the success of our hybrid 
model, as in Ref. [23] we will also explore two other quite different forms for the energy loss 
rate dEjdx, one inspired by perturbative calculations of radiative energy loss and the other by 
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perturbative calculations of collisional energy loss. These expressions are given by 


dE 


o aiJj 

^rad^ ^ 5 

^ collisional 



( 2 . 2 ) 


radiative 


where we shall again treat /^^ad and Ai:coii as parameters to be fixed by fitting to data. This oversim¬ 
plified treatment of radiative and collisional energy loss is not meant to supersede other much more 
sophisticated analysis of these mechanisms [50, 93-98]. Our only goal is to use them as benchmark 
expressions for dEjdx with very different dependence (or lack thereof) on x and to that found 
in (2.1). 

Note that we shall take (2.1) and (2.2) as rates of energy loss for quarks in the parton shower; 
gluons lose more energy. In the case of (2.2), gluons lose more energy by a factor of Ca!Cf = 9/4 
and we shall take /^^ad and FcoW to be larger by this factor for gluons in the parton shower. In the 
case of (2.1), Xgtop is shorter for gluons than quarks by a factor of [64], and as in 

Ref. [23] we shall therefore take /^sc to be larger for gluons in the parton shower than for quarks by 
a factor of (9/4)^/^. 

The energy loss rates Eqs. (2.1) and (2.2) were derived to describe the degradation of the 
energy of partons traversing a static plasma with some constant temperature T. In reality, the 
strongly coupled liquid created in a heavy ion collision is finite in spatial extent and expands and 
cools rapidly, meaning that T depends strongly on both position and time. We will nevertheless 
use these expressions for the instantaneous dE/dx, to describe the energy lost by the partons 
in a shower as they traverse the expanding cooling strongly coupled matter and ultimately the 
energy lost by the jets that emerge from the debris of the heavy ion collision. We describe our 
implementation of this approach in the next Section. 

2.2 The Effects of Flow on the Rate of Energy Loss 

In the previous Section, we have specified the rate of energy loss for an energetic parton traversing 
static plasma with some constant temperature T. In order to use these expressions in the description 
of how an energetic parton loses energy as it traverses the expanding, flowing, cooling, plasma 
created in a heavy ion collision, we will exploit the fact that the expansion, flow and cooling of 
this fluid are described well by nearly inviscid hydrodynamics. This implies that the dynamic 
medium can be described as a collection of fluid cells that are each close to thermal equilibrium, 
locally. As is standard in fluid mechanics, local thermal equilibrium should be understood from a 
coarse-grained point of view: at every fluid cell there is a macroscopic system of size much larger 
than any microscopic scale, such as the inverse temperature of the cell. From this coarse-grained 
perspective, the temperature and the velocity of the fluid cell change from point to point, and from 
time to time. 

Accounting for the variation in the temperature of the fluid is straightforward, and was already 
incorporated in our previous publication [23]. We assume that the temperature T appearing in the 
formulae (2.1) or (2.2) for dEjdx varies in space and time, and at each point in spacetime passed 
by the energetic parton is given by the temperature of the fluid, in the local fluid rest frame, at that 
point in spacetime. The basic assumption behind this adiabatic prescription is that the length scale 
on which an infinitesimal energy loss occurs is small compared to the length scale over which T 
changes. In Ref. [23], we took the variation in temperature into account by integrating (2.1) or 
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(2.2) along the trajectory of each parton in the shower in the collision center-of-mass frame, taking 
T at each point along the trajectory from the hydrodynamic solution describing the bulk fluid. 

The prescription employed in Ref. [23] does not take the velocity of the hydrodynamic fluid 
fully into consideration, as in this prescription the velocity of the moving fluid only affects dE/dx 
in so far as it results in changes in the temperature of the fluid. There is an additional effect, that 
we chose to neglect for simplicity in Ref. [23] but shall incorporate here. It is seen most simply 
by considering a fluid with a constant temperature flowing with a uniform velocity. It is clear 
from their derivations that (2.1) or (2.2) describe the rate of energy loss of an energetic parton 
moving through this fluid in the local fluid rest frame. If this is not the collision center-of-mass 
frame, dE/dx in that frame must be obtained via a Lorentz transformation. In this paper we will 
incorporate this effect into our model description and, in so doing, will improve upon our previous 
treatment of the effects of the medium dynamics on the loss of energy of partons in a jet. Not 
surprisingly, we will find that incorporating the effects of fluid flow on the rate of energy loss has 
significant effects on our results at large rapidity, since it is at large rapidity that the boost between 
the local fluid rest frame and the collision center-of-mass frame becomes large. 

Let us denote the rate of energy loss in the local fluid rest frame by 

^=TFixF,Efj, (2.3) 

where the function E^) is given by the right-hand side of (2.1) or (2.2) and where we 

have highlighted in the notation that dEp/dxp depends upon the distance xp that the parton has 
travelled in the local fluid rest frame and, in the case of (2.1), upon the initial energy E^ that the 
parton had when it was produced at a splitting point in the shower, again as evaluated in the local 
fluid rest frame. In making this statement, we have assumed that the effects of the spatial and 
temporal gradients in the fluid on dEp/dxp can be neglected.^ We now Lorentz transform the 
rate of energy loss (2.3) back to the collision center-of-mass frame, obtaining a result that we shall 
denote by 


(2.4) 

where t is the time in the collision center-of-mass frame since the parton was produced and E[^ is 
the energy that the parton had, in that frame, when it was produced. We can change from x io t 
at will because throughout our treatment we are assuming that the energetic partons in the shower 
move at (very close to) the speed of light. The functions Ef and E are related explicitly by a 
Lorentz transformation that we perform in Appendix A. The result takes on the surprisingly simple 
form 


T{x, Ei^) = Ef{xf, ^in(^))) (2.5) 

^The effects of spatial and temporal gradients in the fluid on the rate of energy loss of an inflnitely heavy quark 
moving through strongly coupled plasma have been computed, to lowest order in fluid gradients [99-101]. They can 
be signiflcant early in a collision, before hydrodynamization. Once the fluid is hydrodynamic, these effects are small. 
We expect that the same is true for the effects of fluid gradients on the rate of energy loss of a massless parton moving 
through the plasma, but we defer checking this by explicit calculation to future work. 
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where and xp are the initial energy and the path length in the local fluid rest frame. These are 
related to quantities defined in the collision center-of-mass frame by 


-Ein = -Ein 7i7' (1 - wv) , (2.6) 

XF{t) = [ dtJ [w2 + (v2 - 2vw + (vw)2)] , (2.7) 

Jto ^ 

where w = p/£’ is the parton velocity, v and are the fluid velocity and Lorentz factor, to 
the time the parton was produced and t is the observation time, all in the collision center-of-mass 
frame. The derivations of these expressions are also given in Appendix A. 

The result (2.5) implies that if the rate of energy loss does not depend explicitly on the energy 
of the parton or the distance that the parton has travelled through the medium, as in the case of the 
collisional rate in (2.2), the fluid velocity will have no effect on dEjdx. In this case, our treatment 
is equivalent to the simpler treatment of Ref. [23]. In the case of the radiative energy loss rate in 
(2.2), or for the strongly coupled rate of energy loss (2.1) that we employ in our hybrid model, we 
expect that including the effects of fluid flow on the energy degradation of jets will be particularly 
important for jets at large rapidity. 

As a simple but illustrative example, let us consider the energy lost by an energetic parton 
propagating through a fluid that is experiencing Bjorken flow, namely boost-invariant longitudinal 
expansion with no transverse flow. If the parton has a large rapidity, the fluid that it is propagating 
through has a large longitudinal velocity meaning that there is a substantial boost between the local 
fluid rest frame and the collision center-of-mass frame. Assuming that both the parton and the boost 
invariant fluid are produced at the same time, the longitudinal velocity of the parton coincides with 
the fluid velocity at its location. If the parton travels a distance L in the collision center-of-mass 
frame during a time f, then according to Eq. (2.7) the distance that the parton travels through the 
fluid in the local fluid rest frame is 

rt 

XF{t)= |wr|rft=—(2.8) 
Jo cosn y 

where and y are the transverse velocity of the parton its rapidity respectively. For particles with 
significant rapidity y, say ^ > 1, we see that the distance they travel in the local fluid rest frame 
is substantially less than the distance they travel in the collision center-of-mass frame. If dEjdx 
grows with x, not taking this Lorentz contraction effect into account, as in Ref. [23], results in an 
overestimate of the amount of energy the parton loses. Conversely, incorporating this effect, as we 
do in this paper, will reduce the energy loss of partons with significant rapidity relative to that in 
Ref. [23]. If the rate of energy loss is given by the expression in (2.2) inspired by radiative energy 
loss, this is the principal effect of flow on the energy loss suffered by an energetic parton. If the 
rate of energy loss is given by the strongly coupled form (2.1), as in our hybrid model, this effect is 
important but it is also important to note that E\^ is also different in different frames. Using (2.6), 


j^p ^ 

cosh^ 


(2.9) 


in this simple Bjorken flow. The effect of this diminution in E^ on the rate of energy loss (2.1) 
is complex. On the one hand, reducing E^ reduces the rate of energy loss. On the other hand, it 
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reduces Xgtop which increases the rate of energy loss. The net effect is not clear a priori. In the 
next Section, we will address this issue by performing a complete simulation of jets in plasma and 
compare it to the results we obtained in our previous publication [23], where we took into account 
the changing temperature along the parton trajectory but left out the effects of the fluid velocity. 

2.3 The Effects of Flow on Single-jet and Dijet observables 

We have re-analyzed the single jet and dijet observables studied in our previous publication [23] 
to include flow effects as prescribed by Eq. (2.5). As already mentioned, the effect of flow for the 
different energy loss models we have studied is different, so we have done this reanalysis for all 
three expressions for the parton energy loss in (2.1) and (2.2). 

Following the same procedure as in Ref. [23], we embedded dijets generated by Pythia 
8.183 [86] into boost-invariant hydrodynamic simulations of the collision dynamics. Initially, as 
in [23], we employed the ideal hydrodynamic simulations of Pb-Pb collisions with center-of-mass 
energy ^/s = 2.76 AGeV computed by Hirano, Huovinen and Nara [89]. These HHN simulations 
incorporate an equation of state obtained from the lattice calculations of Ref. [102]. This equation 
of state has a crossover transition in a range of temperatures 180 < Tc < 200 MeV, somewhat 
higher than the crossover transition temperatures obtained from more recent lattice calculations. 
The value of the crossover temperature Tc has quantitative effects on our results because in our 
model we quench jets by applying an energy loss prescription, (2.1) or (2.2), to the partons in a 
shower only as long as those partons are at points in spacetime where the medium has a temperature 
greater than Tc. When the medium temperature experienced by a given shower parton drops below 
Tc, we turn energy loss off. We make this sharp distinction for simplicity; it is of course unrealistic. 
We expect that using expressions derived for QGP like (2.1) or (2.2) overestimates the energy loss 
at temperatures just above Tc, since the rate of energy loss in a hadron gas is less than that in QGP 
meaning that as the crossover temperature is approached from above the rate of energy loss must 
drop in a way not captured by (2.1) or (2.2). Of course, turning the energy loss off completely is 
certainly an underestimate at temperatures just below Tc. If we turn energy loss off at a specified 
Tc as we do, varying the value of Tc gives us a way to estimate the systematic errors introduced 
into our model results by our crude treatment of late-time energy loss. 

In Ref. [23] we presented all our results with both Tc = 180 MeV and Tc = 200 MeV. We 
do the same in the first two columns of Table 1. The first column repeats results from Ref. [23] 
that include only the effects of variations in T with no effects of the flow velocity on the energy 
loss included. The second column comes from our present calculations, including the effects of 
flow as in Eq. (2.5). In each cell in the table, we quote a range of values for the parameter n^c 
or /^rad or ^coii the expression for the rate of energy loss dE/dx of the partons in the shower 
that we select, (2.1) or (2.2). The range corresponds to fitting Raa for jets with 100 GeV < Px < 
110 GeV to the upper and lower limits of the experimental uncertainty on the measured value of 
this quantity. We then take the entire range spanned by varying Raa for jets with 100 GeV < Px < 
110 GeV and —2<r]<2 over its experimentally allowed range [103] and by varying Tc from 180 
to 200 MeV as the allowed range for a given /^. As described in full in Ref. [23] and reviewed in 
Appendix B, we then compute a variety of single jet and dijet observables over wide ranges in px 
and centrality and obtain good descriptions of the data, in particular using our hybrid model with 
dE jdx given by (2.1). This does not change when we incorporate the effects of flow; what changes 
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Parameter 

HHN hydro 
without flow effects 

HHN hydro 
with flow effects 

SH Hydro 
with flow effects 

Tc r; 

180 MeV 

inge 

200 MeV 

Tcti 

180 MeV 

inge 

200 MeV 

Tc r; 

145 MeV 

fflge 

170 MeV 

^sc 

0.26 - 0.31 

0.30 - 0.35 

0.39-0.46 

0.45-0.53 

0.32 - 0.37 

0.35 - 0.41 

^rad 

0.81 - 1.2 

1.0- 1.6 

1.6-2.4 

2.1 -3.3 

0.97-1.5 

1.2 - 1.8 

^coll 

2.5-3.5 

2.9-4.2 

2.5-3.5 

2.9 -4.2 

1.8-2.6 

2.2-3.0 


Table 1. Values of the fit parameters needed in the specification of dE/dx in our three different energy loss 
models, Eqs. 2.1 and 2.2. The parameters are extracted by comparing model predictions for i?AA for jets 
with 100 GeV< px <110 GeV in central Pb-Pb at ^/s = 2.76 ATeV collisions at the LHC to experimental 
data [103]. The parameters are extracted by employing three different treatments of the bulk hydrodynamic 
fiuid. In the first two columns, we use hydrodynamic simulations from Ref. [89], that we denote by HHN. 
The first column repeats results from Ref. [23]; the second column shows the effects of including the effects 
of the flow velocity on the rate of energy loss as described in Section 2.2. The different choices of Tc 
denote the temperature below which we turn off parton energy loss in our model. In each cell in the table, 
the range of the fit parameter is determined by fitting to one experimentally measured R^a data point, as 
described in the text. In the third column, we use hydrodynamic simulations from Ref. [87] as described in 
the text that were provided to us by Shen and Heinz (SH), with choices of Tc that span the range of crossover 
temperatures favored by current lattice calculations. 


are the fitted values of the parameters, as shown in the first two columns of Table 1. Note that 
^coii does not change since, as we discussed in the previous Section, flow does not affect dEjdx 
if dEjdx is given by the collisional expression from (2.2). Both /^sc and increase: in these 
cases, incorporating the effects of flow reduces the energy loss of partons with substantial rapidity 
and since we are fitting to Raa for jets with pseudorapidities over the range — 2 < p < 2 we have 
to increase and A^rad- 

Throughout the rest of this paper, we will embed jets and implement our hybrid model for 
their quenching in a new set of hydrodynamic simulations of Pb-Pb collisions at ^/s = 2.76 ATeV 
and 5.02 ATeV provided to us by Shen and Heinz, based upon the codes developed by them and 
their collaborators in Refs. [87] and [88]. These 2-1-1-dimensional simulations, with boost-invariant 
longitudinal expansion, incorporate an equation of state from Ref. [104] (referred to there as the 
s95p-vl-PCE150 equation of state) that incorporates results from lattice QCD calculations and 
from a hadron resonance gas at low temperatures. With this equation of state, it is possible for us 
to set the Tc below which we stop parton energy loss within a range 145 MeV < Tc < 170 MeV 
that reflects results for the range of crossover temperatures from current lattice calculations, see for 
example Refs. [105, 106]. For simplicity, the hydrodynanfic calculations employ a temperature- 
independent Tj/s — l/(47r). Also for simplicity, we employed ensemble-averaged MC-Glauber 
initial conditions [107], neglecting all effects of event-by-event fluctuations. Because we are not 
concerned with the evolution below Tc, we evolved the initial conditions using viscous hydrody¬ 
namics alone, with no cascade afterburner. 

As shown in the third column of Table 1, we have refitted the parameters a^sc, ^rad and a^coII 
in the expressions (2.1) and (2.2) to the experimentally measured Raa for jets with —2<r]<2 
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and 100 GeV< pj <110 GeV using the SH hydro simulations and setting Tc to 145 MeV or 170 
MeV to yield some sense of the theoretical systematic uncertainty as before. The values of the 
parameters so obtained are smaller than those in the second column, reflecting the lower value of 
Tc meaning that energy loss continues for the partons in the showers for a longer time. We have 
checked that most of the difference between the parameter values in the second and third columns 
is indeed due to the change in the value of Tc, as anticipated. Perhaps coincidentally, the parameter 
values in the third column are quite similar to those from Ref. [23] given in the first column: the 
consequences of adding flow effects and lowering Tc to a large degree cancel. 

In Appendix B, we update the results of Ref. [23] for the jet Raa as a function of p^ and 
centrality, and for the dijet imbalance in Pb-Pb collisions at ^/s = 2.76 ATeV, comparing them to 
CMS data [6, 103]. We see that with SH hydro and a lower Tc, just as with HHN hydro and a higher 
Tc in Ref. [23], when we fit the one parameter n^c in our hybrid model to one data point the model 
does a very good job of describing these full data sets. In Appendix C, we present results for Raa 
and the dijet imbalance for Pb-Pb collisions with ^/s = 5.02 ATeV. These constitute predictions 
for the LHC heavy ion run coming late this year. 

The improvements we have introduced in our description of the dynamics of the expanding 
droplet of QGP in going from the first column of Table 1 to the second column to the third column 
have resulted in changes to the numerical values of /^sc. ^rad and However, the conclusions 
about the implications of the values of these parameters that we reached in Ref. [23] all remain 
unchanged. The extracted values of /^rad and a^coII imply such large “weak” couplings or such large 
logarithmic corrections to the leading order energy loss result that a perturbative analysis is called 
into question.^ The extracted value of /^sc m the strongly coupled expression (2.1) for dEjdx in 
our hybrid model corresponds to a stopping distance for energetic partons in the strongly coupled 
QGP of QCD that is about three to four times longer than that in the strongly coupled plasma 
of A/" = 4 SYM theory. And, as can be seen by comparing the plots in Appendix B to those in 
Ref. [23], the hybrid model with the improvements that we have introduced continues to provide 
a good description of data on both single jet and dijet observables over a wide range of jet p^ and 
event centrality. 

2.4 Species Dependence of Jet Suppression 

Although the central conclusions we draw from confronting our hybrid model with the measured 
observables compiled in Appendix B are unchanged from those we reached in Ref. [23], the conclu¬ 
sions that we drew there about the species dependence of jet suppression are substantially affected 
by the inclusion of the effects of flow on parton energy loss that we have described in the preceding 
Sections. In the left panel of Fig. 1, we repeat a calculation from Ref. [23], now with SH hydro and 
the lower Tc range from the third column of Table 1, and reach different conclusions. We show the 
ratio of the number of quark-initiated jets to the number of gluon-initiated jets with jet momenta in 
the range 100 < Px < 150 GeV and \r]\ < 2 as a function of centrality for the three models we 
have studied. As in [23], we have defined the species of the jet-initiator from our Pythia “data” 
(in a way that is impossible to do in experimental data) as the identity of the parton that initiated 
the DGLAP shower to which the hardest particle in the jet belongs. The inclusion of the effects of 

^For an analysis of the effects of large logarithmic corrections to medium parameters, see Refs. [108, 109]. 
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Figure 1. Left panel: Ratio of quark-initiated jets to gluon-initiated jets with \r]\ < 2 as a function of 
centrality for different models. Right panel: Ratio of quark-initiated jets to gluon-initiated jets as a function 
of jet pseudorapidity for p-p collisions in vacuum, according to Pythia. 


flow on parton energy loss has resulted in the bands corresponding to our three different models 
being much closer to each other in Fig. 1 than they were in Ref. [23]. The power to discriminate 
between models of energy loss implicit in Fig. 1 is much less than was the case in Ref. [23]. This 
change arises as a combination of two facts. First, as we discussed in Section 2.2 the inclusion of 
the effects of flow reduces the energy loss of partons (both quarks and gluons) at large rapidity if 
dE/dx is given by the strong coupling expression (2.1) or the radiative expression from (2.2), but 
leaves things unchanged if the collisional expression from (2.2) is employed. Second, as shown in 
the right panel of Fig. 1, without any energy loss there are many more quark-initiated jets at large 
rapidity than at mid-rapidity for jets with the same range of as in the left panel of the Figure. As 
the pseudorapidity increases, jets in this interval of transverse momentum become more and more 
quark dominated. The consequences of these two facts in concert is that the many quark-initiated 
jets at large rapidity lose less energy once flow effects are included than they did in the analysis of 
Ref. [23], if dE/dx depends on flow as in either the strongly coupled or radiative expressions. In 
these cases, the effects of flow on the rate of energy loss serve to enhance the relative weight of 
the large rapidity jets in the final jet sample, and these large rapidity jets are predominantly quark- 
initiated. This means that, after summing over \r]\ < 2 as in the left panel of Fig. 1, the inclusion of 
the effects of flow pushes the strongly coupled and radiative Nq/Ng bands in this Figure upwards, 
bringing both much closer to the collisional band. This happens even though, as we mentioned 
in Section 2.1 and discussed in detail in Ref. [23], the Casimir-dependence of the rate of energy 
loss in (2.1) is weaker than in (2.2): this difference in the Casimir dependence of dE/dx pushes 
the strongly coupled Nq/Ng band in the left panel of Fig. 1 downward relative to the collisional 
band, but the effects of flow in concert with the rapidity-dependence of the production of quark 
and gluon jets in p-p collisions brings it back up. Consequently, our current results render the 
species-dependence of quenching rather insensitive to the underlying energy loss dynamics. 

We close this section by noting that despite the substantial rapidity dependence of the species- 
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Figure 2. Pseudorapidity dependence of jet Raa in our hybrid model for ^/s = 2.76 ATeV (left) and 
= 5.02 ATeV (right). 


dependence of jet production, and of jet quenching, the resulting net jet Raa is remarkably inde¬ 
pendent of rapidity in our hybrid model, as shown in Fig. 2. This was noted previously in Ref. [110], 
and we find the same behavior in our control models in which we choose the dependence of dE jdx 
according to either one of the expressions (2.2). Although given the strong species-dependence 
manifest in the right panel of Fig. 1 the flatness of the rapidity dependence seen in Fig. 2 appears 
coincidental, it is in agreement with experimental results from the ATLAS collaboration [17]. 


3 Boson-Jet Correlations, Including Predictions for ^/s = 5.02 ATeV Collisions and 
for Z-jet Correlations 

The extensive exploration in Ref. [23] demonstrated that the hybrid model that we have developed 
describes currently available inclusive jet and dijet data from LHC heavy ion collisions rather 
successfully. In Appendix B, we compile updated versions of these results, using SH hydro and 
taking flow effects into account as we described in the previous Section, and we reconfirm the 
conclusions of Ref. [23]. This gives us considerable confidence in the model framework. However, 
this class of observables proves not to be very sensitive to whether we choose dEjdx as in (2.1) 
or (2.2) [23]. The distinctions between the strongly coupled form for dE/dx and our two control 
models provided by these observables are not sufficient to differentiate between these different 
hypotheses for the microscopic dynamics of energy loss, at least with present uncertainties. The 
strongly coupled form for dEjdx does provide a better description of the dijet imbalance, but 
its predictions are not sufficiently distinct from those of the control models. We must, therefore, 
consider further observables. 

In deciding how to go further, we face a choice. We could start adding more physics to the 
model, which would allow us to confront new classes of observables. For example, if we were 
to add in the transverse momentum kicks that the medium delivers to the shower partons passing 
through it, and the recoil that the shower partons delivers to the medium, we could engage with jet 
shape observables. This is an attractive prospect, but we defer it to future work. Before adding to 
the model, and in particular before adding a second free parameter, it is our responsibility to first 
ask whether there are further observables that our present model, with its single parameter already 


- 14 - 




fixed as described in the previous Section, could reasonably be expected to describe. This is our 
goal in this paper. 

For most of the remainder of this paper, we turn to observables involving jets back-to-back 
with either an energetic and isolated photon or a Z-boson. We will keep the parameter /^sc in our 
hybrid model fixed exactly as in the previous Section, and confront our model with existing data on 
three different observables describing 7 -jet correlations in Pb-Pb collisions at ^/s = 2.76 ATeV. As 
before, our model describes the data well, although it is fair to note in advance that the statistical 
error bars on these measurements are substantial because to this point the data sets of 7 -jet events 
have not been large. Nevertheless, the confrontation of the hybrid model with these three new 
sets of data provides a strong independent validation of our hybrid model, which has no further 
adjustable parameters. 

Our central purpose in this Section is to use our hybrid model to provide predictions for the 
measurements that the data sets with higher statistics by about an order of magnitude that are 
anticipated late this year will make possible. To this end, we provide predictions for three 7 -jet 
observables and three Z-jet observables in Pb-Pb collisions at = 5.02 ATeV. In making these 
predictions, we will keep the parameter /^sc in the hybrid model and the parameters /^rad and a^coII 
in the control models set to the same values that we obtained in the previous Section by fitting 
to one data point at y/s = 2.76 ATeV. The principal change in going from y/s = 2.76 ATeV to 
y/s = 5.02 ATeV is a (modest) increase in the temperature of the plasma. The principal effects of 
this increase arise from the explicit temperature dependence in ( 2 . 1 ) or ( 2 . 2 ), and so are included 
in our analysis. We are leaving out any small decrease in the values of the /^’s that may arise if the 
plasma becomes slightly less strongly coupled, anticipating that this effect will likely be too small 
to be resolved given other uncertainties, both theoretical and experimental. 

Although the observables that we study in this Section are similar in some respects to those 
that describe dijet events, since neither 7 ’s nor Z-bosons interact strongly with the plasma, when 
one triggers on an energetic 7 or Z-boson the jets produced in the same hard scattering have points 
of origin that are distributed through the collision volume differently than is the case for dijet pairs, 
where selection bias favors points of origin closer to the surface where at least one of the jets will 
suffer less energy loss. This means that the jets produced in association with a 7 or a Z-boson 
sample a different path-length distribution than the jets in dijet pairs. A second difference between 
boson-jet events and dijet events is that both 7 ’s and Z-bosons are much more likely to be produced 
in association with a quark jet than with a gluon jet. The third difference is that since the 7 or Z- 
boson suffers no energy loss, as they do not interact with the medium, their energy and direction are 
good proxies for the initial energy and initial direction of the hard parton (usually, the quark) going 
in the opposite direction. This is quite different than in the case of dijet events, where generically 
both jets should be expected to have lost some energy. Of course, the hard parton opposite a 7 or 
Z-boson may split into more than one jet; in the analysis of 7 -jet and Z-jet correlations, it cannot 
be assumed that the boson is back-to-back with only one jet in the final state. 

3.1 Generation and Selection of Monte Carlo Events 

For our study of photon-jet correlations we analyzed 10^ hard scattering processes, in Pythia 
8.183 [86] p-p collisions at ^/s = 2.76 TeV and another 10^ such events in Pythia p-p collisions 
at = 5.02 TeV. We require at least one photon with a transverse momentum above a desired 
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cut, typically choosing = 60 GeV, and with a pseudorapidity in the range \r]j\ < 1.44. 

We set the parameter in Pythia (basically, the minimum momentum transfer in the hard 
processes that Pythia is sampling) safely lower than for example choosing px™^^ = 40 GeV 
if we are recording events that include a photon with > 60 GeV, since Pythia does not reliably 
reproduce the photon spectrum all the way down to its We placed each of these 10^ events 

at random in the transverse plane of a Pb-Pb collision, with ^/s = 2.76 ATeV or 5.02 ATeV, 
choosing the location according to the probability distribution for the number of binary collisions 
occurring at each transverse position according to an optical Glauber model. We keep the direction 
of the photon (in azimuth and in rapidity) as it was in the Pythia p-p collision from which we 
have taken it. We generate the hard processes in Pythia without underlying p-p events because 
we embed each hard process into a hydrodynamic description of the matter produced in the heavy 
ion collision, meaning that including the underlying p-p event from Pythia would be double¬ 
counting. As we described in the previous Section, we use the boost-invariant 2+1-dimensional 
hydrodynamical simulations of the expanding cooling droplet of matter produced in the heavy ion 
collision from Ref. [87] to determine the local temperature and local velocity vector of the fluid at 
the position of each of the partons in the jet shower described by Pythia, as they propagate through 
the cooling hydrodynamic fluid as a function of time. We follow each parton in the shower, and as 
long as a parton sees a local temperature that is greater than Tc we reduce its energy at a rate dE/dx 
given by (2.1) in our hybrid model, or by (2.2) if we are investigating one of our two control models. 
When a parton that has lost some energy branches, we start each of its daughters off with the same 
fractional energy loss as the parent had when it branched. The daughters then lose further energy 
according to (2.1) in our hybrid model or (2.2) in our control models. In this way, we compute 
the energy lost by each of the partons in the showers produced by the hard scattering process, in 
so doing modifying the partons in the final state relative to what they were in the original Pythia 
event. As in Ref. [23], we neglect energy loss in the pre-equilibrium stage of the collision before 
the hydrodynamic evolution of Ref. [87] starts at r = 0.6 fm. We have described these procedures 
in greater detail in Ref. [23]. As in the previous Section, we vary our choice of the temperature 
Tc below which we turn off the energy loss over the range 145 < Tc < 170 MeV. The result of 
our analysis is 10^ hard scattering events (for each collision energy and for each choice of Tc) in 
which the partons accompanying the photon have lost energy as a result of their passage through 
the hydrodynamic medium. 

Next, we analyze the 10^ events in order to “measure” each observable that is of interest to 
us. In doing so, we only analyze the (modified) parton showers coming from the hard scattering. 
We do not include any of the particles that would be created as the hydrodynamic fluid freezes out. 
This means that we do not have any background particles and we therefore do not perform any 
background subtraction in our calculation of observables. When we compare to data, we of course 
compare to measurements made after the experimentalists have done their background subtraction. 
For this reason, for those observables which have not been fully unfolded in the experimental 
analysis we include a smearing procedure, to be describe in more detail later, to mimic resolution 
effects. 


"^For one of our observables, we wanted to include photons with px^ > = 40 GeV in our analysis. For this 

observable, we ran 10^ events withpx’^'’^ = 20 GeV. 
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Although we of course know (from Pythia) exactly which photons in our events are the pho¬ 
tons coming from the originating hard process (prompt photons) and which photons are instead 
produced in the parton shower (fragmentation photons), we do not use this information in our com¬ 
putation of observables. Instead, we perform an isolation analysis of all high energy photons in our 
Monte Carlo data, patterned upon what experimentalists must do. This procedure allows us to study 
a sample of photons which is closer to that in experimental analyses, since in experiments prompt 
and fragmentation photons can only be distinguished via isolation cuts, with finite efficiencies. As 
it turns out, quenched jets associated with prompt photons constitute about 80% of our final “data” 
sets, with about 20% coming from fragmentation of jets in events without a prompt photon. 

For readers who would like a little more detail, our procedure is as follows. We select events 
containing at least one photon with and \r]^\ < 1.44. After quenching the parton 

showers in these events as described above, following the analysis of Refs. [8, 111] we consider a 
photon to be isolated if the sum of the energy of all the particles within a cone around this photon 
of radius Ar = y^A0^~+~A^ = 0.4 is below 5 GeV. Contrary to what is done in the analysis 
of experimental data, we perform this isolation cut at the partonic level, instead of at the hadronic 
level. In the unlikely situation that two or more isolated photons are found in a single event, we 
treat the one with the highest transverse momentum as the leading photon. Next, we construct 
a sample of photon-jet pairs from our sample of isolated photons. We reconstruct jets using the 
anti-ZcT algorithm [112, 113], typically choosing R = 0.3, and keep those events in which we find 
an associated jet at an angular distance A(/) > Ttt/S from the isolated photon. This angular cut 
suppresses the contribution of events in which the isolated photon is a fragmentation photon and 
events with more than one associated jet to our final photon-jet sample. 

We generate a sample of Z-jet events via an analogous procedure. We select 10^ Pythia 
events in which a hard scattering process produces a Z-boson with p^^ > again typically 

choosing = 60 GeV. We need not worry about fragmentation Z-bosons: the large mass of 
the Z makes it very unlikely that they are produced in the parton showers. We therefore need 
not apply any isolation procedure: any Z is a prompt Z. For the sake of simplicity, we use the 
same p^ and p cuts applied to photon events. We reconstruct jets using the anti-/cT algorithm as 
before, again typically choosing i? = 0.3, and obtain our Z-jet sample by requiring an associated 
jet at an angular distance Acj) > Ttt/S from the Z. In experimental measurements, Z-bosons are 
reconstructed through their di-muon decay. Since high energy muons are not modified at all by 
their passage through the plasma, and since the properties of the associated jets are the same no 
matter how the Z-boson decays, to increase the size of our sample we simply do not allow our 
Z-bosons to decay, keeping all of them in our sample. 

3.2 Photon-Jet Observables: Comparison with Experimental Results at = 2.76 ATeV 
and Predictions for = 5.02 ATeV 

We use the Monte Carlo samples of photon-jet events, prepared and quenched as described above, 
to construct different observables which can be confronted with experimental data. In this section, 
we describe these observables, we show the results obtained from our model for collisions at ^/s = 
2.76 ATeV and compare them to data, and we provide the predictions that we obtain from our 
model for collisions at ^/s = 5.02 ATeV. 
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Figure 3. Distribution of the transverse momentum imbalance of photon-jet pairs, xj^ = for 

Pb-Pb collisions. The left and right panels show the 0-10% and 10-30% most central events, respectively. 
The upper panels show our results for collisions with y/s = 2.76 ATeV, as well as data from Ref. [111]. The 
lower panels show our predictions for heavy ion collisions with ^/s = 5.02 ATeV, anticipated for late this 
year. 


We first construct the photon-jet imbalance. This is analogous to the dijet imbalance, except 
that here we have a photon instead of a leading jet. In the case of dijets we define the associated 
jet as the jet in the pair that has less transverse momentum, meaning that the ratio of the transverse 
momentum of the associated jet to that of the leading jet is less than one by definition. This is not 
so in the case of the photon-jet imbalance, since the associated jet can have more or less transverse 
momentum than the photon. Defining xj^ = , this observable can be less than or greater 

than one. For example, in p-p collisions in vacuum, an associated jet can have more transverse 
momentum than the photon if there is a second jet in the event, in the same hemisphere as the 
photon. In a heavy ion collision, the energy loss experienced by the partons propagating through 
the medium pushes xj^ downwards. However, if the passing jet sweeps particles from the plasma 
into the jet cone this can in principle push xj^ upwards, but this effect is expected to be small 
at large energies. Our hybrid model neglects this possibility, meaning that in every event in our 
sample xj^ is less than (or equal to) what it would have been in the absence of the medium. 

In Fig. 3, we show the distribution of the imbalance in the transverse momentum of the as¬ 
sociated jet relative to that of the photon, xj^, for events with two different centralities and for 
events with two different collision energies. Following the conventions established in the experi¬ 
mental analyses in Refs. [8, 111], these distributions are normalized to the total number of photons, 
rather than to the total number of photon-jet pairs; integrating each of the curves in Fig. 3 therefore 
yields a number below one. Also as in the experimental analysis, we only consider photons with 
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transverse momentum > 60 GeV. The associated jet is reconstructed with the anti-ZcT algo¬ 
rithm with R = 0.3; we only count events in which the associated jet has a transverse momentum 
> 30 GeV. The widths of the colored bands that illustrate our results in this Figure — and in 
many Figures that follow — incorporate both the uncertainty that comes from varying Tc between 
145 and 170 MeV and the uncertainty that comes from varying our model parameter /^sc over its 
allowed range 0.32 < /^sc < 0.41, determined in Section 2. That is, we obtain four curves by 
repeating our calculation of the observable in question, here xj^, with Tc and /^sc each set to its 
lowest and its highest value, and plot the band that extends from whichever one of the four curves 
is lowest to whichever curve is highest at each point in the Figure. 

As a reference, in Fig. 3 we display the xj^ distribution in p-p collisions at the same nucleon- 
nucleon energy as predicted by Pythia, i.e. with no medium-induced parton energy loss. In order 
to mimic the effects of jet-energy resolution on the transverse momentum of the associated jet ob¬ 
tained in the analysis of the Pb-Pb data after the subtraction of the background, in our p-p results we 
have smeared the momenta of the associated jets obtained from our Monte Carlo calculation with a 
centrality-dependent Gaussian broadening. The parameters of this smearing, reported in Ref. [114], 
were tuned to reproduce the measured distributions after embedding PYTHlA-generated photon-jet 
samples into real lead-lead events. For this reason, the smeared proton-proton distribution is not 
identical in the left and right panels of Fig. 3, which display our results for Pb-Pb collisions with 
two different centralities. For the present, these smeared p-p results are the correct reference to 
which both the experimental results as measured by the CMS collaboration [7] and our results 
for quenched jets in Pb-Pb collisions should be compared. (See Ref. [115] for fully unfolded ex¬ 
perimental results to come.) We apply the same smearing procedure to our simulated quenched 
jets, represented by the colored band. The smearing parameters are at present known only for 

= 2.76 ATeV, not for ^/s = 5.02 ATeV. For the present, we have decided to employ the 
same smearing parameters at the higher collision energy. Our predictions therefore assume that the 
effects of background subtraction on these observables are similar at these two collision energies. 

The differences between the (smeared) proton-proton Monte-Carlo data and the results of 
our in-medium calculations displayed in Fig. 3 are due to energy loss. The sample of isolated 
photons we have used to construct these distributions is dominated by prompt photons, which do 
not lose momentum when traversing the plasma.^ The partons in the showers that become the 
associated jets, however, interact strongly with the medium and lose energy according to (2.1) as 
they propagate through it. In some cases, this pushes the transverse energy of the associated jet 
below 30 GeV, meaning that the event does not get counted as a photon-jet in our Pb-Pb analysis 
although it was counted in our p-p analysis. This is why the integrals of the curves illustrating our 
Pb-Pb results in Fig. 3 are smaller than the integrals of our p-p results. More on this below. In 
other cases, when the transverse energy of the associated jet remains above 30 GeV, the effect of 
the energy loss is to reduce displacing the photon-jet imbalance distribution toward smaller 
values of xj^. 

Keeping in mind that we fixed the single parameter in our hybrid model by comparing it to 
the single-jet suppression Raa^ see Section 2, it is remarkable how well the photon-jet imbalance 

^The small fraction of fragmentation photons which fulfill the isolation requirement do suffer energy loss via the 
quenching of the partons from which they originate. This effect is small for this observable, but it can have consequences 
for more differential observables, as we will discuss below 
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Figure 4. Fraction of events with an isolated photon in which we find a photon-jet pair, which is to say in 
which we find an associated jet with > 30 GeV at an azimuthal angle more than Ttt/S away from that 
of the isolated photon. We plot this fraction as a function of the photon transverse momentum, in collisions 
with ^/s = 2.76 ATeV (left) and ^/s = 5.02 ATeV (right). The colored band shows the results from our 
hybrid model, with its strongly coupled form for the rate of energy loss. For comparison, the violet dots 
show the smeared p-p calculations (see text for details). Our results for y/s = 2.76 ATeV are compared to 
CMS data [111]. 


distribution that we obtain from our model agrees with CMS measurements in ^/s = 2.76 ATeV 
collisions in both centrality bins in Fig. 3. The fact that only one side of photon-jet pair loses 
energy makes the interpretation of this observable cleaner than in the case of the dijet imbalance. 
Of course, at present the large statistical uncertainties in the photon-jet measurements is a limitation 
on their use to differentiate between different model assumptions for the rate of energy loss dEjdx, 
which is to say a limitation on their use as diagnostics of the mechanism of energy loss. We 
illustrate this point in Appendix D, where we present the results that we obtain by repeating our 
photon-jet and Z-jet calculations using the control models for dEjdx in (2.2). We find that these 
control models make predictions that are distinct from those of our hybrid strong/weak coupling 
model, with dEjdx given by (2.1), but the distinctions are too small to be resolved by the present 
data, with their statistical uncertainties. 

Given these considerations, perhaps the most important results that we can use our model 
to provide are our predictions for the upcoming runs of the LHC, at ^/s — 5.02 ATeV, where 
the photon-jet data sets are expected to be larger by about an order of magnitude and the statis¬ 
tical uncertainties are expected to be substantially smaller than at present. For this reason, in the 
lower panels of Fig. 3 we show the predictions of our hybrid model, with its strongly coupled 
form for the rate of energy loss, for the photon-jet imbalance distribution in Pb-Pb collisions at 
^ = 5.02 ATeV. We have applied the same kinematical cuts used in current measurements (and 
our calculations) at ^/s = 2.76 ATeV. The imbalance distribution shows little dependence on the 
collision energy. 

To further test the success of our hybrid approach, we now turn to exploring other photon-jet 
observables. In Fig. 4 we show the fraction of isolated photons that come with an associated jet, as 
reconstructed with the anti-/cT algorithm with R = 0.3, that has 30 GeV and A(/) > TnjS. 
We plot this quantity, which we denote by Rj^, as a function of the transverse momentum of the 
isolated photon for heavy ion collisions with ^/s = 2.76 ATeV (left) and ^/s = 5.02 ATeV (right). 
In the plots in Fig. 4, the smeared proton-proton Pythia simulations are represented by the violet 
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Figure 5. Ratio of the transverse momentum spectra of jets associated with an isolated photon in Pb-Pb 
collisions to that in p-p collisions, as a function of the jet transverse momentum, for two different ranges 
of the photon transverse momentum and for = 2.76 ATeV (upper panels) and = 5.02 ATeV (lower 
panels). The hybrid model with its strongly coupled form for the rate of energy loss describes the available 
CMS data at ^/^ = 2.76 ATeV well. 


dots, the results of our hybrid strong/weak coupling model are represented by the colored band, 
and the experimental results from Ref. [Ill] for ^/s = 2.76 ATeV are the black data points.^ The 
broad xj^ distribution seen in Fig. 3 implies that this ratio must be an increasing function of for 
both collision energies, since the 30 GeV cut on the associated jet energy is more and more easily 
satisfied as the momentum of the photon is made larger and larger. Quenching reduces Rj^ since 
it pushes the energy of some of the associated jets below 30 GeV. As for the photon-jet imbalance 
distribution, we find good agreement between the Rj^ obtained from our hybrid model and that 
measured in present experiments with y/s = 2.76 ATeV, in this case for all values of the photon 
momenta.^ Our predictions for heavy ion collisions with = 5.02 ATeV are shown in the right 
panel of Fig. 4. As before, this observable shows little sensitivity to the collision energy, at least 
within these kinematical cuts. 

We have also analyzed the spectrum of jets produced in association with an isolated pho¬ 
ton. In Fig. 5, we show the ratio of the spectrum of associated jets in Pb-Pb collisions to that in 
proton-proton collisions, Iaa, for two different ranges of photon transverse momenta, 60 GeV< 

^Note that the ratio of the integral of the photon-jet imbalance in Pb-Pb collisions to that in p-p collisions, see Fig. 3, 
is the ratio of Rj^ for all photons with pR > 60 GeV in Pb-Pb collisions to that in p-p collisions. The fact that the 
colored band lies below the violet dots in Fig. 4 was therefore foreshadowed in Fig. 3. 

^Note that in order to extend our calculations of Rj^ down to 40 GeV< pR <50 GeV we used a sample of Monte 
Carlo events with the Pythia parameter set to 20 GeV. 
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<80 GeV and > 80 GeV. The observable Iaa^ and in particular its suppression below 
1, can be thought of as the photon-jet analogue of the single-jet suppression Raa, but instead of 
being constructed for inclusive jets as Raa is, Iaa is constructed using only the associated jets 
back-to-back with an isolated photon. This implies that the distribution of the jet energies as well 
as the fragmentation pattern and hence the jet masses, of the p-p jets and of the Pb-Pb jets that 
enter into the calculations of Iaa and Raa are different. Furthermore, the distribution of the point 
in the transverse plane at which the hard scattering event that produces a jet selected in an Iaa 
analysis is quite different from that for the jets selected in an Raa analysis. It is therefore striking 
that even though we fitted the single parameter in our hybrid model to a single measured value of 
Raa, when we compare the results for the photon-jet Iaa obtained from the model with data at 

= 2.76 ATeV, displayed in the upper panels of In Fig. 5, we see such good agreement, for both 
ranges of the photon energy and over the whole range of jet p^. In the lower panels of Fig. 5, we 
show the predictions from our hybrid model with its strongly coupled rate of energy loss for heavy 
ion collisions with ,/s = 5.02 ATeV. As for the previous photon-jet observables, the spectrum of 
the associated jets hardly changes between these two collision energies. 

The agreement between the predictions of our hybrid strong/weak coupling model, with its 
Pythia branching and its strongly coupled form for the rate of energy loss, for all of these photon- 
jet observables and the data available today is very encouraging. Having fixed the single parameter 
in our model using a single measurement of Raa for inclusive jets, without introducing any new 
parameters we have obtained a good description of the experimental data for a total of 5 different 
observables, one involving inclusive jets, one involving dijets, and three involving photon-jets, all 
with their centrality and energy dependence. These observables sample different in-medium path 
length distributions of the quenched jets, different shapes of the original jet spectrum, different 
fragmentation patterns and jet mass distributions, and different quark vs. gluon compositions of 
the observed jets. Despite all these differences, our model is able to describe the systematics 
observed in all the data correctly. To avoid over-reaching in drawing conclusions, however, it is 
important to explore the predictions of the different control models described in Section 2.1 for 
these observables. The results of this analysis can be found in Appendix D. There are distinctions 
between the predictions of the control models and our hybrid model, but these distinctions are small 
compared to the statistical uncertainties in present data. It is therefore not possible at present to use 
the agreement between our hybrid model and photon-jet data to argue that the data favors a strongly 
coupled form for the rate of energy loss. We therefore await the higher statistics data expected later 
this year in collisions with ^/s = 5.02 ATeV, and have provided the predictions of our hybrid model 
for all three photon-jet observables for collisions at this higher energy. We have also provided such 
predictions for Z-jet observables below, in Section 3.3. In addition, it is important to investigate 
the predictions of our hybrid model for intra-jet observables like fragmentation functions, as we 
shall do in Section 4. These are all paths toward using jet data in the service of understanding the 
dynamics of the interaction between energetic partons and strongly coupled plasma. 

3.3 Z-Jet Observables: Predictions for = 5.02 ATeV 

In this section, we turn to the Z^-jet observables that are complementary to the photon-jet observ¬ 
ables of the previous Section. The large Z-mass ensures that Z-bosons seen in a heavy ion collision 
were almost without exception produced promptly, in hard processes dominated by short-distance 
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Xjzo 

Figure 6. Distribution of the transverse momentum imbalance of Z-jet pairs, xjz = for Pb- 

Pb collisions with ^/s = 5.02 ATeV. The left and right panels show the 0-10% and 10-30% most central 
events, respectively. Here and below, the colored bands show the results from our hybrid model, with its 
strongly coupled form for the rate of energy loss and the violet dots show the smeared p-p calculations for 
comparison. 


physics: it is extremely unlikely for a Z-boson to be produced during the fragmentation of a parton 
in a jet, unless the jet energy is much higher than is relevant to us. The Z-boson production mech¬ 
anism is therefore under good theoretical control.^ Furthermore, because of their large width and 
short lifetime, Z-bosons decay during the very early stage of a heavy ion collision, even prior to 
plasma formation. The Z-bosons that are identified in heavy ion collisions are those that decayed 
leptonically, in particular via This means that their decay products do not inter¬ 

act strongly with the pre-equilibrium partonic matter or with the strongly coupled plasma, once it 
forms. Z-bosons are in this respect similar to prompt photons, making Z-jet events similar in their 
utility to photon-jet events. The Z-jet events come with the added advantage almost all Z’s are 
prompt Z’s. As consequence, the experimental identification of Z-bosons via their decays 

do not require isolation cuts, which leads to an arguably cleaner determination of the associated jet 
energy. The only disadvantage of Z-jet events is that Z’s are less numerous than photons. 

We construct the same class of observables for Z-jet events that we constructed for photon- 
jet events in the previous Section. Because of the low statistics of Z-boson production in y/s = 
2.76 ATeV collisions, no constraining measurements exist at present for the Z-jet observables that 
we shall construct.^ We therefore present predictions for these observables in the ^/s = 5.02 ATeV 
heavy ion collisions coming soon in LHC heavy ion Run 2. We have chosen the same kinemat- 
ical cuts for the Z-jet observables that we (and the CMS collaboration) have used in photon-jet 
observables, to facilitate comparison between our results for the two cases. Of course, once the 
experimentalists decide on the cuts that they will use for their Z-jet analyses of the data-to-come, 
we can re-run our analyses with their cuts. 

In Fig. 6, we show the distribution of the Z-jet imbalance observable xjz = in 

heavy ion collisions with two different ranges of centrality with ^/s = 5.02 ATeV In both cases the 
colored band shows the predictions of our hybrid model, with its strongly coupled form (2.1) for 
the rate of energy loss, for Pb-Pb collisions and the violet dots show the distribution of xjz for p-p 

^Note, of course, that modifications of the nuclear parton distribution function in the relevant intermediate a;-region 
can alter the production rate of Z-bosons, by modifying the composition of the initial fiux of partons. 

^For preliminary low-statistics measurements of Z-jet correlations at = 2.76 ATeV, see Ref. [116]. 
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Figure 7. Fraction of events with a Z-boson in which we find a Z-jet pair, which is to say in which we find 
an associated jet with > 30 GeV at an azimuthal angle more than 77r/8 away from that of the Z-boson, 
in collisions with ^/s = 5.02 ATeV. 


collisions as predicted by Pythia. As we did in our analysis of photon-jet observables, we have 
smeared the momenta of the associated jets in both our Pb-Pb and p-p calculations. For the present, 
before better guidance becomes available once experimentalists have begun the analysis of LHC 
Run 2 data, we have used the same smearing functions here as we (and the CMS collaboration) 
used for the photon-jet observables that we discussed in the previous Section. We obtained the xjz 
distributions in Fig. 6 from a sample of events in which we required a Z-boson with > 60 GeV 
and an associated jet reconstructed using the anii-kx algorithm with R = 0.3 that has > 
30 GeV and is separated in azimuthal angle from the Z-boson by Acf) > Ttt/S. We required that 
both the Z and the associated jet have I 77 I < 1.6. As in the case of the xj^ distribution in photon-jet 
events, the xjz distribution is broad in p-p collisions, indicating the importance of events with a 
Z and two jets, in particular those arising from initial state radiation. This means that at present 
Z-bosons are not substantially better as taggers of the associated jet energy than photons are, which 
motivates the future development of methods to suppress events with more than one jet in the final 
state. Because the transverse momentum of the Z-boson and its decay products do not 

change in the medium, the difference between the Pb-Pb distribution and the p-p distribution in 
our calculation is entirely due to the energy lost by the partons in the associated jet in the Pb-Pb 
collisions due to their passage through the strongly coupled plasma. As for photon-jet events, we 
see a reduction in the integral of xjz, analyzed further below, and a displacement of the distribution 
toward smaller xjz- The magnitude of this displacement is comparable to the corresponding shift 
in the xj^ distribution in photon-jet events, see Fig. 3. As in that case, there are clear but small 
distinctions between the results we obtain with our hybrid model, shown in Fig. 6 , and those we 
obtain when we use our control models for the rate of energy loss (2.2) instead. We present these 
in Appendix D. 

In Fig. 7, we compute the fraction of Z-bosons in our sample that come with an associated 
jet, as reconstructed with the anii-kx algorithm with R = 0.3, that has > 30 GeV and 
A(j) > Ttt/S. We plot this quantity, which we denote by Rjz, as a function of the transverse 
momentum of the Z-boson. We show the results obtained from our smeared proton-proton PYTHIA 
simulations as the violet dots and the predictions for Pb-Pb collisions from our hybrid strong/weak 
coupling model as the colored band. As in the photon-jet case, as the partons in the associated jet 
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Figure 8. Ratio of the transverse momentum spectra of jets associated with a Z-boson in Pb-Pb collisions 
with ^/s = 5.02 ATeV to that in p-p collisions as a function of the jet transverse momentum, for two different 
ranges of the Z-boson transverse momentum. 

shower lose energy the total energy of the associated jet can drop below our = 30 GeV cut, 
meaning that energy loss leads to a reduction of the Z-jet yield as a fraction of the number of Z- 
bosons. This makes the integral under the colored bands in Fig. 6 less than that under the violet dots 
there, and it pushes the colored band in Fig. 7 below the violet dots. The qualitative behavior and 
magnitude of this reduction is comparable to the reduction in the analogous photon-jet observable, 
see Fig. 4. 

Finally, in Fig. 8 we analyze the spectrum of jets produced in association with a Z-boson. We 
show the ratio of the spectrum of associated jets in Pb-Pb collisions to that in p-p collisions, Iaa^ for 
two different ranges of Z-boson transverse momenta, 60 GeV< <80 GeV and >80 GeV. 
The effect of energy loss on this observable is again comparable to its effects on the analogous 
photon-jet observable, see Fig. 5. 

In Appendix D, we repeat the analysis of Rjz and Iaa for jets produced in association with a 
Z-boson in our control models, where we use the expressions (2.2) for the rate of energy loss. We 
find that these observables exhibit little sensitivity to the form of dE/dx, meaning little sensitivity 
to the microscopic dynamics via which the partons in the jet shower interact with the strongly 
coupled plasma. 

In Sections 3.2 and 3.3, we have provided predictions for three photon-jet observables and 
three Z-jet observables in heavy ion collisions with ^/s = 5.02 ATeV obtained from our hybrid 
strong/weak coupling model, with fragmentation taken from Pythia and a rate of energy loss 
(2.1) as at strong coupling. We await the data that will come from LHC heavy ion Run 2 with 
considerable anticipation. 

4 Fragmentation Functions 

We now turn to the analysis of a more differential class of jet observables, namely fragmentation 
functions. We saw in Ref. [23] and have confirmed in Fig. 17 of Appendix B that the predictions of 
our hybrid model and of our two control models for the ratio of the partonic fragmentation function 
of inclusive jets in PbPb collisions to that in p-p collisions are distinct for the three models. This 
motivates the hope that the higher statistics measurements expected from the coming LHC run 
may serve to distinguish between models. And, it motivates us to compute the predictions of 
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all three models for the partonic fragmentation functions of jets produced in association with an 
isolated photon or a Z-boson. We shall present the results of these calculations in Section 4.1. 
We were initially surprised to see that the partonic fragmentation function ratio for jets produced in 
association with bosons turn out to be more similar for the three models we are considering than was 
the case for inclusive jets. Understanding this effect, which we shall do by the end of Section 4.1, 
leads us, in Section 4.2, to introduce a new observable constructed from the fragmentation functions 
of jets in dijet pairs in Pb-Pb collisions, without the need for any p-p reference. This turns out to be 
the observable that is most effective at differentiating among our three models of any observable 
that we have considered to date. 

Unlike the more inclusive observables described in the previous Section, the hadronic frag¬ 
mentation functions that experimentalists measure are quite sensitive to hadronization effects. 
However, the dynamics of hadronization, even in vacuum, are not under full theoretical control. 
In fact, the predictions of different Monte Carlo event generators for fragmentation functions 
can differ among themselves, and in comparison with p-p data, by as much as 20% as shown 
in Refs. [117, 118]. Since the modification of the fragmentation functions in Pb-Pb collisions 
with respect to those in p-p collisions is itself on the order of several tens of percent at most (see 
the measurements reported in Ref. [9, 15, 119]) and since there are differences in hadronization 
dynamics (in particular, differences in the patterns of color recombination) for jets in Pb-Pb and 
p-p collisions [120-122], it will be challenging to compare the ratios of partonic fragmentation 
functions in Pb-Pb collisions to those in p-p collisions — ratios which we calculate in our models 
in Section 4.1 — to data. The in-medium effects of interest are comparable in magnitude to the 
known uncertainties coming from our lack of understanding of hadronization dynamics. This is 
strong further motivation for the importance of the observable that we introduce in Section 4.2: the 
ratio of the fragmentation function of inclusive jets in Pb-Pb collisions to that of the associated jets 
in dijet pairs. Since this is the ratio of the fragmentation functions of two different classes of jets in 
Pb-Pb collisions, with no need for a p-p reference, many hadronization uncertainties will cancel. 

We shall restrict our calculations to partonic fragmentations throughout this Section, meaning 
that our calculations are not sensitive to hadronization and so are not affected by its challenges. 
The observable that we introduce in Section 4.2 is the one in this Section for which this will to 
the greatest degree possible also be true in experimental data. For this observable, as for those we 
calculate in Section 4.1, however, the fragmentation functions describing fragments with the lowest 
Pt’s cannot be described reliably by our hybrid model or by the two control models because none 
of these models include the contribution to the low-px component of a jet arising from the wake in 
the recoiling plasma that the jet plowing through it produces [47, 123-127]. 

4.1 Fragmentation Functions of the Associated Jets in Photon-Jet and Z-jet Pairs 

Fragmentation functions are defined as the distribution of hadrons within a jet with a given fraction 
z of the total longitudinal momentum py of the jet. Longitudinal, here, means in the direction of the 
jet axis. In p-p collisions, where most of the activity in events with hard jets comes from hadrons 
produced via the fragmentation and subsequent hadronization of the virtual partons produced in 
an initial elementary partonic collision, fragmentation functions provide us with information about 
how the showering process via which the large virtuality of the initial partons relaxes takes place. 
Furthermore, since final state effects are negligible, py obtained from the energy of a jet recon- 
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structed with a sufficiently large reconstruction radius provides a good proxy of the initial energy 
of the hard parton that fragmented to form the jet. This means that the 2 ;-fraction of the final frag¬ 
ments are directly related to one of the QCD evolution variables. This is not the case in Pb-Pb 
collisions. 

We have already mentioned one of the complications in Pb-Pb collisions: at low some of 
the particles in a reconstructed jet did not originate from the initial hard parton that was produced 
in the initial elementary partonic collision. Some soft particles reconstructed in the jet come, in¬ 
stead, from the hadronization of moving quark-gluon plasma, set in motion by the momentum that 
the jet passing through it transfers to the medium through which it is passing. Operationally, this 
enters the analysis because in analyzing jets in Pb-Pb collisions it is necessary and standard to do 
a background subtraction to remove the hadrons formed from the quark-gluon plasma, and this 
background subtraction procedure is based upon the assumption that the momenta of these hadrons 
is uncorrelated with the direction of the jet. To the extent that this is the case, the background 
subtraction removes, on a statistical basis, particles in the jet cone that are not part of the jet itself. 
However, since the interaction of the jet with the plasma transfers momentum to the plasma, this 
back reaction (or recoil) effect means that the background subtraction procedure cannot remove 
all the particles from the medium: there is no way to disentangle all of them from the products of 
the jet shower; some of them must end up incorporated into jet observables. This means that, in 
Pb-Pb collisions, fragmentation functions at low are sensitive to physical processes other than 
jet fragmentation. Addressing these additional physical processes requires a dynamical treatment 
of the response of the medium to the passage of the jet, which is beyond our current model imple¬ 
mentation. For this reason, the results of our calculations become less reliable at small z. In this 
section, we will look at jets with p^ > 30 GeV, meaning that py > 30 GeV. So, ln(l/ 2 :) = 2.7 or 
ln(l/ 2 :) = 3.5 corresponds to fragments with p^ > 2 GeV or p^ > 1 GeV. We will plot our results 
out to larger values of ln(l/ 2 ;), smaller values of z, but the effects of medium recoil that we are not 
including become more and more important for ln(l/ 2 :) > 3. In Section 4.2 where we consider 
inclusive jets and dijets, rather than jets produced in association with a boson, we will look at jets 
with Pt > 80 GeV, meaning that our results there will be reliable out to somewhat larger values of 
ln(l/z). 

There is a second complication in interpreting jet fragmentation functions in Pb-Pb collisions, 
and this is that the py of a jet, as reconstructed from the final state, is less than the energy of the 
initial hard parton that fragmented into the jet because the partons in the jet have lost energy as 
they propagate through the strongly coupled medium produced in the Pb-Pb collision. At least 
some of the “lost” energy (according to the data [6, 128], a significant amount of it) ends up as soft 
particles moving at large, random, angles relative to the jet axis, and is not included when the jet 
is reconstructed. Therefore, even when they are reconstructed with a large reconstruction radius, 
the total energy of the quenched in-medium jets is smaller than the total energy of the hard partons 
originating from an elementary partonic collision. This means that when a jet is reconstructed in a 
Pb-Pb collision, it is impossible to make an experimental determination of the energy of the initial 
hard parton. Consequently, if one constructs a fragmentation function using the standard definition 
of z, namely the ratio between the momentum of the momentum of an individual hadron to the 
py of the whole jet as reconstructed, this 2 ; is not directly related to the evolution variable in a 
DGLAP shower and, more generally, is simply hard to interpret. We will nevertheless report the 
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predictions of our hybrid model and the two control models for fragmentation functions computed 
in this standard way, for comparison to future data. 

Because we are looking at jets produced in association with a photon or Z-boson, however, 
there is an obvious alternative. We can define Zr^ or (we shall denote these variables generically 
as 2 ;^) as the ratio of the momentum of an individual hadron in the jet to the momentum of the 7 or 
Z-boson, using the momentum of the electroweak boson in the event (which cannot have lost any 
energy since it does not interact with the quark-gluon plasma) as a proxy for the momentum of the 
initial hard parton that later fragmented and lost energy, forming the jet that the experimentalists 
reconstruct^^. This is an improvement but it is not a panacea: we have seen in the previous Section 
that even in p-p collisions there is a broad distribution of the momentum imbalance between a 
boson (photon or Z) and the associated jet in a boson-jet event. One significant contributor to 
this imbalance is the fact that in many events the elementary hard scattering process produces one 
photon or Z plus more than one hard parton, not just one. Regardless, this imbalance is entirely due 
to perturbative vacuum QCD physics, not to any in-medium effects. This means that even when 
we construct fragmentation functions using or we are not reliably dividing by the actual 
momentum of the initial hard parton that fragmented into the jet we are looking at. Still, by using 
z^ or 2 ;^ in Pb-Pb collisions we are using a variable that is as good a proxy for what we want as is 
the case in the standard fragmentation function in p-p collisions. 

We first compute the fragmentation functions constructed with respect to the reconstructed 
energy of the jet produced in association with a photon or Z-boson. (We use the standard variable 
2 :, but here we denote it zj to emphasize that its denominator is the momentum of the reconstructed 
jet. It is defined by zj = In Fig. 9 we plot the ratio of the fragmentation functions 

in Pb-Pb to p-p as obtained in our framework for jets of px > 30 GeV produced in association 
with isolated photons (upper panels) and Z-bosons (lower panels) with px > 60 GeV. We have 
performed these simulations for two different Pb-Pb centralities, 0 —10 %, shown on the left panels 
and 10 — 30 %, shown on the right. Both photons and Z-bosons are required to have I 77 I < 1.4 while 
the jets are constrained to |77jet| < 1.6. The energy and direction of the jet axis are reconstructed 
with the anti-ZcT algorithm with radius i? = 0.3. Following the experimental analyses of inclusive 
fragmentation functions in Ref. [15], these fragmentation functions are constructed by including 
all particles surrounding the jet axis within an angular distance in ij-cj) space of r = 0.3. Since 
our simulations do not include the underlying event, these correspond to particles from the hard 
scattering process that remain correlated with the jet direction. (As noted above, in the analysis 
of experimental data a background subtraction is done, but soft particles from the plasma that is 
in motion following the jet will not be subtracted. This means that our calculations should not 
reproduce the data-to-come at large values of ln{l/zj).) 

The general features of the fragmentation function ratios plotted in Fig. 9 are very similar to 
those that we found for inclusive jets in Ref. [23] and have confirmed in Fig. 17 of Appendix B. 
All the models display an enhancement of the hardest part of the fragmentation function in Pb-Pb 
collisions relative to p-p collisions. This is a generic behavior of any mechanism that removes 
soft particles from the jet, either via energy loss as here or via deflecting them into a direction 
far from that of the jet [129]. Removing soft particles increases the fraction of jets with a few 

See Ref. [47] for a similar definition of the scaling variable. 
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Figure 9. Predictions of our hybrid model, with strongly coupled energy loss for the partons in a Pythia 
shower, as well as our two control models for the partonic fragmentation function ratios (fragmentation 
function for jets in Pb-Pb collisions over that for jets in p-p collisions) for jets produced in association with 
an isolated photon (upper panels) or a Z-boson (lower panels) in Pb-Pb collisions at a/s = 5.02 ATeV at two 
centralities (left and right panels). The fragmentation functions are constructed with respect to the variable 


zj 


parton / Jet 

P|| /Fj| . 


hard fragments, which leads to the increase in the hard part of the fragmentation function. For 
all models there is also a depletion in the Pb-Pb fragmentation function at intermediate zj. This 
is the expected result from quenching, which tends to reduce the energy of the fragments that 
propagate in plasma. Remarkably, for the energy range of bosons and jets explored in those figures, 
and with our current uncertainties, the pattern of fragmentation at large and intermediate zj is 
indistinguishable among the three models we explore, despite their very different path length and 
energy dependences. We comment further on this below. At smaller values of zj, the hybrid model 
with its strongly coupled energy loss suppresses soft fragments more than the control models. 
However, this separation between models occurs in a regime where the fragments have momenta 
smaller than 2 GeV, meaning that our calculations of fragmentation functions are not reliable there. 
Adding in the contributions from a medium that has picked up momentum from the jet passing 
through it, meaning that it is not completely removed by the background subtraction, would push 
the Pb-Pb fragmentation functions up in this soft region by an amount that our model does not 
permit us to estimate at present. 

In Fig. 10 we reanalyze the fragmentation function ratios, this time using the boson momentum 
Zb to define the scaling variable according to zb = — with the transverse 


- 29 - 





























Figure 10. Predictions of our hybrid model and our two control models for the partonic fragmen¬ 
tation function ratios (Pb-Pb over p-p) for jets produced in association with a boson as a function of 
Zb = — for two different centralities for photon-jet (upper) and Z-jet (lower) events 

in Pb-Pb collisions at y/s = 5.02 ATeV. 


momentum of the isolated photon or Z-boson.^^ As already mentioned, our main motivation for 
redefining the scaling variable is to have a better proxy for the jet energy prior to quenching. If 
all bosons were prompt, their momentum would be insensitive to in-medium effects and energy 
loss would only affect the numerator of zb> For those prompt bosons, the mismatch between the 
boson and the initiator parton of the jet originates entirely in vacuum processes. This is the case for 
Z-jet correlations, where the possibility of producing a Z in a jet shower is highly suppressed. In 
contrast, in the isolated photon sample we use, there is a small fraction of fragmentation photons 
even after we make our isolation cut, and for these photons energy loss effects are present via the 
quenching of their parent parton. Therefore, for photons, there is also a small dependence of the 
denominator of on quenching; the effect of this small dependence on the fragmentation function 
ratio depends on the zb range. 

Note that as a consequence of the broad distribution of the boson-jet imbalance even in p-p 
collisions, zb can be larger than one: there are jets with larger than the momentum of the boson 
they are associated with. Those events populate the region of negative \ogl /zb in Fig. 10. In 
this region, all models lead to an enhancement of the in-medium fragmentation function of jets 

Because of the fluctuations of the rapidity of the centre of mass of the elementary partonic collision that leads to 
boson-jet events, the rapidity of the boson and the jets do not need to be correlated, unlike the transverse momentum. 
For this reason, to construct zs we have chosen to use only transverse momenta and chosen not to project the momenta 
of fragments along the boson direction. 
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associated to photons. At first this seems puzzling, but there are in fact two reasons for this effect. 
The first is the fact that jet quenching does reduce the energy of fragmentation photons, as their 
parent partons lose energy. The second reason comes from imposing an isolation cut on photons 
in the events in our calculations, which do not include the particles corresponding to the medium. 
Since quenching affects all partons in our events, prompt photons in events with more than one jet 
are more isolated in our Pb-Pb simulations than in vacuum. This leads to an enhancement of the 
fragmentation function ratio for jets produced in association with isolated photons in the region of 
xj > 1. However, it is not clear to us whether this effect will persist in a full simulation of photon- 
jet events in which particles from the Pb-Pb background are incorporated in the sample. We leave 
the study of the fate of this enhancement for future work. For the Z-jet correlation neither of these 
two effects are present, and the negative log zb region is slightly suppressed, as expected. 

The ratio of fragmentation functions also exhibits a non-trivial structure in the vicinity of 
log 1/zB 0.5 {zb 0.6). This structure is correlated with the position of the maximum of the 
in-medium boson-jet imbalance in Figs. 3 and 6, xj^ ^ 0.6 and xjz ^ 0.6. Since the maximum 
of the in-medium imbalance distribution is shifted towards smaller x values, the non-monotonic 
behavior of the fragmentation function ratio in this region reflects the behavior of the imbalance 
distribution. Indeed, the relative abundance of associated jets with xj^ ^ 0.6 or xjz ^ 0.6 is 
enhanced in Pb-Pb collisions with respect to p-p. The fact that all three models exhibit the same 
behavior is a consequence of the coincidence of the imbalance distribution in the three models. 

As for fragmentation functions constructed with the variable zj, here too the strongly coupled 
model only separates from the control models for soft particles where the physics of how the 
medium responds to the passage of the jet, physics that none of our models includes, becomes 
important. 

In summary, the analysis of the fragmentation functions of jets produced in association with 
photons and Z-bosons indicates that mechanisms of energy loss that do not increase the number of 
hard fragments in jets, like the ones we have explored, lead to robust modifications to the fragmen¬ 
tation pattern of these jets. This also means that this type of observable is not very sensitive to the 
microscopic mechanism of parton energy loss. 

We close this section by recalling our initial motivation — the separation between the pre¬ 
dictions of our models for the fragmentation functions of inclusive jets seen in Ref. [23] and in 
Fig. 17 — and asking why that separation between model predictions is less in the fragmentation 
functions for jets produced in association with photons and Z-bosons seen in Figs. 9 and 10. The 
answer comes in understanding the selection effects in a sample of “inclusive jets”. Because the 
jet production spectrum is a steeply falling function of pj, and because usually two or more jets 
are produced in an event, most of the jets in a sample of inclusive jets are the most energetic jet in 
an event. In selecting a sample of inclusive jets, one is preferentially selecting jets that are the jet 
in their event that has lost the least energy. In contrast, when one selects jets by first identifying 
an isolated photon or Z-boson and then reconstructing an associated jet there is no such selection 
effect. This means that, on average, the inclusive jets whose fragmentation functions are shown in 
Fig. 17 have travelled through the medium over a shorter path-length and, again on average, they 
are jets that fragmented less. Fewer fragments, i.e. a jet with a lower jet mass and a smaller open¬ 
ing angle and a harder fragmentation function, means that within the jet there are fewer partons 
losing energy in the medium, and therefore means less energy loss. Both these effects are likely 
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small compared to the event-by-event variation. But, on average, inclusive jets contain somewhat 
fewer fragments^^ and traverse somewhat less plasma. What are the consequences in our models 
of the fact that in going from Fig. 17 to Figs. 9 and 10 there is an increase in path length? The 
path-length dependence of dE/dx in (2.1) is stronger than in the control models (2.2) meaning that 
in the mid-range of 2 ; where we see the effects of quenching (say 1 < log 1/z < 3) we expect the 
increase in path length to push the predictions of our hybrid strongly coupled model down relative 
to the control models. What are the consequences of the fact that in going from Fig. 17 to Figs. 9 
and 10 the jets become wider and there is some reduction in the energy of the partons within them? 
In (2.1), a reduction in the energy of the partons means a reduction in the stopping length Xgtop 
and an increase in dE/dx. So, for this reason also we expect to see the predictions of our hybrid 
strongly coupled model pushed down relative to those of the control models. And, this is indeed 
what we see. Unfortunately, the predictions of the hybrid model get pushed down just to the extent 
that in Figs. 9 and 10 they are essentially on top of the predictions of the control models. 

In the next Section, we show that it is possible to select a sample of jets in which the pre¬ 
dictions of the hybrid strongly coupled model for the fragmentation functions in the mid-range 
1 < log 1 /z < 3 are pushed down even farther. 

4.2 Fragmentation Functions of the Associated Jets in Dijet Pairs 

Motivated by the results and discussion above, we now turn our attention to the fragmentation 
functions of the associated (less energetic) jets in dijet events. Whereas in the inclusive jet sample 
of Fig. 17 we have selected jets that are likely to be those among the jets in their event that have 
been quenched the least, by selecting associated jets in dijet events we will likely be selecting those 
that have been quenched the most. That means we will be selecting those that have, on average, 
traversed a longer path-length of medium and those that were produced, on average, with a larger 
jet mass and jet opening angle and that, on average, contained more, and therefore lower energy, 
fragments before quenching. 

In the left panel of Fig. 11, we show the ratio of fragmentation functions for jets of 80 < 
Pt < 300 GeV with 0.3 < I 77 I < 2 produced in association with a jet of > 120 GeV in Pb- 
Pb collisions to that in p-p collisions. By the associated jet in a dijet pair we will always mean 
the jet with the lower energy. For comparison, in the right panel of Fig. 11 we have analyzed 
inclusive jet fragmentation in a lower energy range, such that the momentum of the inclusive jets is 
in the same range as the momentum of the associated jets in the left panel of Fig. 11. Although by 
careful comparison to Fig. 17 we see that the modification pattern seems to be slightly dependent 
on the jet energy, the basic features in the right panel of Fig. 11 and in particular the ordering of 
the predictions of the three models, is the same as in Fig. 17. The inclusive jets are selected in a 
way that makes them likely to be the less quenched jets in their event with, on average, a shorter 
path length and fewer fragments. The jets in Section 4.1, selected via having been produced in 
association with photons or Z-bosons, have no such selection effects. And, the associated jets in 

^^There are actually two reasons why they contain somewhat fewer fragments. First, as we shall see only in the next 
Section, even in p-p collisions the higher energy jet in a dijet pair — and this is what most jets in an inclusive jet sample 
are — tends to have a harder fragmentation function. And, second, jets with fewer, harder, fragments lose less energy 
[130] in the plasma produced in a Pb-Pb collision and so are even more likely to end up being the higher energy jet in a 
dijet pair. 
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Figure 11. Left: Model predictions for the Pb-Pb to p-p ratio of fragmentation functions of associated jets 
(lower energy jets) in dijet events at ^/s = 5.02 ATeV with a leading jet of px > 120 GeV. Right: Same 
for inclusive jets with the same range of px and 77 as the associated jets in the left panel. The right panel is 
similar to Fig. 17, but is for a different px range and is for the higher LHC heavy ion Run 2 collision energy. 
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Figure 12. Model predictions for the ratio of the fragmentation functions of inclusive jets to the frag¬ 
mentation functions of associated jets in dijet events whose most energetic jet has px > 120 GeV in Pb-Pb 
collisions at a/s = 5.02 ATeV. Both the inclusive jet and associated jet samples are constrained to the same 
77 and Px range, see the text. We see a striking, and robust, separation between the predictions of our hybrid 
strongly coupled model and our two control models. We also show the same ratio, constructed with the same 
kinematic cuts, for inclusive and associated jets in p-p collisions. The ratio in p-p collisions should not be 
compared to that in Pb-Pb collisions: each stands on its own and has its own implications. (Comparisons 
between Pb-Pb and p-p are shown in Fig. 11.) Here, the p-p results illustrate the differences between jets se¬ 
lected inclusively and jets selected by association with a leading jet in vacuum, with no jet quenching. And, 
the comparison between the Pb-Pb results from the different models shows how the interplay between these 
differences and the path-length and energy dependence of different expressions for dEjdx yields predictions 
for this ratio that depend sensitively on the underlying microscopic dynamics of jet quenching. 


the left panel of Fig. 11 are selected in a way that makes them likely to be the more quenched jets 
in their event with, on average, a longer path length and more fragments. Sure enough, we see that 
the predictions of our hybrid strongly coupled model — with a dEjdx that depends strongly on 
path-length and that increases at lower energies as Xgtop decreases — are pushed lower than those 
of the control models in the left panel of Fig. 11. 

The model-dependence of the fragmentation function ratios seen by comparing the left and 
right panels of Fig. 11 is striking: the ordering of the predictions of the three models is opposite 
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in the two panels, with the hybrid strongly coupled model predicting the least depression of the 
fragmentation function of the inclusive jets and the most depression of the fragmentation function 
of the associated jets. Seeing this motivates us to introduce a new observable in Fig. 12 in which 
we, in effect, take the ratio of the right panel of Fig. 11 to the left panel of Fig. 11. That is, 
we propose to compare the fragmentation function of inclusive (which is to say leading) jets to 
the fragmentation function of associated (lower energy) jets in Pb-Pb dijet events. To avoid trivial 
kinematic differences between the energies of the inclusive and associated jets, we have constructed 
this ratio with jets in the same interval of and rapidity, 80 < px < 300 GeV and 0.3 < \r]\ < 2. 
As in Fig. 11, the associated jets whose fragmentation functions constitute the denominator in the 
ratio plotted in Fig. 12 are produced in association with leading jets of > 120 GeV. (We have not 
investigated other choices of leading jet p^.) The predictions of our hybrid strongly coupled model 
and of our two control models for this new observable are displayed in Fig. 12. This observable 
yields the largest separation between the predictions of our models of any observable that we have 
investigated. This means that it is particularly sensitive to the underlying microscopic dynamics 
behind the modification of jets in medium. 

The new observable that we have introduced in Fig. 12 has the added virtue that its construc¬ 
tion does not require fragmentation functions from jets in p-p collisions. It is the ratio of the 
fragmentation functions of differently selected Pb-Pb jets with the same kinematics. This means 
that none of the uncertainties coming from the differences in the way hadronization occurs in Pb- 
Pb collisions and p-p collisions come into play. Our calculations behind Fig. 12 are, like those in 
Section 4.1 and Appendix B, calculations of partonic fragmentation functions. Here, though, we 
expect that the ratio of fragmentation functions displayed in Fig. 12 should be a better predictor for 
the same ratio of the hadronic fragmentation functions that experimentalists can measure than in 
the case of any of the other fragmentation function ratios that we have constructed because much of 
the uncertainty in our understanding of hadronization should cancel in this ratio. Of course, in the 
softest region, at the smallest z’s, it continues to be the case that we are leaving out contributions to 
the fragmentation functions from the moving quark-gluon plasma, set in motion by the jet passing 
through it. This effect should at least partially cancel in the ratio plotted in Fig. 12; this is yet one 
more advantage of using the ratio of fragmentation functions of differently selected jets with the 
same kinematics in the same Pb-Pb collisions, without introducing jets from p-p collisions. 

In Fig. 12 we have also plotted (as purple dots) the same ratio of the fragmentation function 
of inclusive jets to that of associated jets for p-p collisions from Pythia, with the same kinematic 
cuts as in our Pb-Pb analysis. We show this solely to confirm that even in the absence of any 
medium there is a difference between the inclusive jets and the associated jets. Selecting inclusive 
jets selects jets with fewer and therefore harder fragments. Selecting associated jets selects jets 
with more and therefore softer fragments. In Pb-Pb collisions, then, in the presence of a medium 
the differences between the rates of energy loss come into play as we have discussed. This effect, 
plus the difference between the path lengths seen by inclusive and associate jets on average, are 
amplified by the different path-length and energy dependence of the rates of energy loss in our hy¬ 
brid strongly coupled model and our two control models. The result is the large separation between 
the predictions of the models for the ratio plotted in Fig. 12, making this ratio so discriminating. 
We look forward to seeing experimental measurements of this ratio; they have great potential to 
teach us about the microscopic dynamics that results in the observed modification of jets in heavy 
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ion collisions. 


5 Conclusions and Outlook 

In this paper we have explored a broad range of jet observables in Pb-Pb collisions at LHC energies 
within the context of the hybrid strong/weak coupling model that we introduced in Ref. [23]. This 
is a phenomenological approach to the physics of jet quenching in which we aim to synthesize the 
very different types of physical processes involved in the production, branching, and subsequent 
in-medium dynamics of jets in heavy ion collisions. In particular, our model separates the short 
distance physics which controls the production and hard branching evolution of jets, which behaves 
as in vacuum, from the longer distance processes that control the interaction of the jet shower with 
the strongly coupled fluid produced in energetic heavy ion collisions. Our goal is to describe the 
short distance physics with standard tools for describing the weakly coupled physics of jets in vac¬ 
uum, in particular with Pythia, and to describe the longer distance processes using insights from 
calculations done using strong coupling methods. In particular, to model the long distance physics 
we adopt results obtained by analyzing the rate of energy loss dEjdx of an energetic massless 
parton propagating through the hot liquid plasma of strongly coupled A/^ = 4 SYM theory, as ob¬ 
tained via gauge/gravity duality. We apply these results in QCD upon assuming that the differences 
between the hot liquid plasma phases of the two theories can be absorbed in a single parameter 
which controls the stopping distance of energetic excitations in plasma. We fitted this one model 
parameter to jet Raa data in our previous publication. Our model, now fully specified, yields a 
very good description of the suite of inclusive jet observables and dijet observables discussed in 
our previous publication and, as our explicit analysis in this paper shows, of all currently available 
data on photon-jet correlations. All these observables in sum span a wide range of energies, orig¬ 
inate from different hard production processes with different primordial spectra, and correspond 
to jet samples with different selection biases, different ratios of quark jets to gluon jets, different 
fragmentation patterns and jet mass distributions, and that traverse different distributions of path 
lengths. Nevertheless, we have obtained a satisfactory description of all of these observables us¬ 
ing our simple, hybrid, implementation of strongly coupled dynamics — namely by taking jets 
from Pythia and applying the strongly coupled rate of energy loss dEjdx parton-by-parton to the 
partons in a parton shower. 

Although at one level we are pleased by the now increasingly many successes of our hybrid 
model approach, at another level they are frustrating. Our model is simple, gluing together two 
rather different descriptions of physics at different scales and in so doing incorporating much un¬ 
known physics into just one parameter. Of course the first goal in creating such a model is to 
capture some of the physics correctly, and it seems that we have done that. However, it could be 
even more interesting to see the model breaking down and to use ways in which it fails to describe 
some feature of some experimental data to understand which of the aspects of the physics that the 
model leaves out are important, and how, and why, and where. The string of successes in the com¬ 
parison between the results of our model calculations and experimental data preclude investigations 
of this nature at present. 

All that said, our explicit analysis of two control models with parametrically different expres¬ 
sions for dEjdx cautions us that we should not rush to conclude from the successes of our hybrid 
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model that the experimental measurements of the various jet, dijet and photon-jet observables fa¬ 
vor the strongly coupled form for dE/dx over other possibilities. Rather, we must acknowledge 
that these observables are not strongly sensitive to the parametric form of the rate of energy loss 
dE jdx. The experience that we have gained from our analysis suggests that as long as the vacuum¬ 
like branching processes that are at the core of jet dynamics are described well, any mechanism that 
is able to quench particles, in particular the softer partons in a jet, can capture the bulk features of 
the measured distributions of jet observables, including the various dijet and photon-jet asymme¬ 
tries and correlations, as long as the parameter that governs the overall magnitude of dEjdx is 
chosen appropriately. 

Nevertheless, it is certainly heartening and perhaps even remarkable that the range of values 
of the one parameter in our hybrid strong/weak coupling model, that we find provides a good 
description of so much jet data agrees so well with a priori expectations, should be smaller 
than but of order one, exactly as we have obtained. A stopping length for energetic partons in the 
strongly coupled QGP of QCD that is three to four times longer than that in the strongly coupled 
J\f = 4 SYM plasma with the same temperature, as we find, is an eminently reasonable result. 

These conclusions are all conclusions that we reached in our previous publication. The analy¬ 
ses of photon-jet observables presented here serves to reinforce them in many ways. 

Looking ahead, we hope that among the many calculations from our hybrid model that we have 
presented the ones that will be most important will be the many predictions that we are making for 
experimental measurements that are anticipated in the near future. These can be grouped into three 
categories, the first two of which are: 

• We have provided the predictions of our hybrid model for inclusive jet observables, dijet 
observables, and photon-jet observables in heavy ion collisions with collision energy ^/s = 
5.02 ATeV, in anticipation of LHC heavy ion Run 2. This is particularly important for the 
photon-jet observables. The present data have low statistics, and correspondingly large error 
bars, making the fact that our hybrid model describes them well less impressive than it could 
be. In the run to come, the statistics will be greater by about an order of magnitude, meaning 
that the error bars should be significantly smaller, making the confrontation between the 
predictions of our model and these measurements much more constraining. 

• We have provided the predictions of our hybrid model for Z-jet observables, which we hope 
will be tested in future LHC runs. 

Confrontation between these predictions and the data to come should serve either to further strengthen 
our confidence in the approach to jet quenching that we have introduced or to identify and quan¬ 
tify ways in which it fails, ideally pointing toward which aspects of the physics that we have left 
out are most important and guiding the improvement of the model. It is also possible that as the 
experimental uncertainties shrink measurements of these observables could serve to differentiate 
between the different assumptions about the dynamics of parton energy loss, and consequent dif¬ 
ferent forms for dE/dx. However, the distinctions between the predictions of our hybrid model 
and our two control models for these observables are small, limiting the discriminating power of 
these observables even as the experimental uncertainties shrink. 

It seems clear that in order to find observables that provide more discrimination among the 
different possible dynamical processes via which the partons in a shower lose energy as they tra- 
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verse the strongly coupled plasma we will need to investigate intrajet observables. Utilizing those 
observables that involve the angular shape of jets must wait, as it will require adding further physics 
to the model, including for example the transverse momentum picked up by the shower partons as 
they interact with the medium, and so must involve the addition of at least one new parameter to 
the model. There is good motivation for such investigations, but we leave them to future work. 
With the one parameter model that we have constructed here, the class of intrajet observables that 
we may be able to describe is those constructed from fragmentation functions: 

• The third category of predictions that we have made are predictions for various ratios of 
fragmentation functions. We have taken advantage of having a tool with which we can com¬ 
pute multiple observables (here, fragmentation functions for jets in p-p collisions, inclusive 
jets and jets in dijet pairs in Pb-Pb collisions, and jets produced in association with pho¬ 
tons or Z-bosons in either Pb-Pb or p-p collisions) in varying kinematic regimes to search 
for new discriminating observables which are particularly sensitive to differences between 
mechanisms of energy loss. We have provided the results of our hybrid model calculations 
for ratios between partonic fragmentation functions in photon-jet and Z-jet events in Pb-Pb 
collisions to those in p-p collisions, but for reasons that we have understood qualitatively 
these turn out not to be particularly discriminating. Also, these predictions are sensitive to 
differences between hadronization dynamics in Pb-Pb and p-p collisions that are not under 
good theoretical control at present, reducing the reliability of our predictions for these ratios. 
The most discriminating observable that we have found is the ratio of the fragmentation func¬ 
tion of inclusive jets in Pb-Pb collisions to the fragmentation functions of jets in the same 
kinematic regime, in the same collisions, that are the lower energy jets in a dijet pair. (We 
refer to the latter as associated jets.) Unlike most jet observables, this measurement does not 
require any p-p reference data. And, since the ratio we propose is constructed from jets in 
the same kinematic regime in the same PbPb collisions, differing only in how they were se¬ 
lected, we expect that many of the theoretical uncertainties associated with the modification 
of hadronization in medium should cancel. This makes the ratio of partonic fragmentation 
functions that we compute a better proxy to the ratio of fragmentation functions that experi¬ 
mentalists will measure. Furthermore, we find that the predictions of our hybrid strong/weak 
coupling model and our two control models for this ratio, see Fig. 12, are well separated over 
a wide range of z, making the discriminating power of this observable robust even after the 
softening of the fragmentation functions expected after hadronization. And, perhaps best of 
all, we have a good qualitative understanding of why this ratio is such a discriminating ob¬ 
servable. First, in Pb-Pb collisions as in vacuum inclusive jets tend to contain fewer, harder, 
fragments than associated jets on average. Second, the distribution of the path length of the 
medium through which a sample of inclusive jets has propagated is, on average, shorter than 
that for a sample of associated jets. Both these effects mean that, on average, the inclusive 
jets have lost less energy than the associated jets. And, we have shown that both these ef¬ 
fects push the predictions of our hybrid strong/weak coupling model and our control models 
apart. Measurement of this ratio of fragmentation functions should be a particularly effective 
way to gain information about the dynamics via which energetic partons lose energy as they 
traverse strongly coupled plasma. The experimental determination of this ratio in the immi- 
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nent LHC heavy ion Run 2 can therefore shed light on the microscopic dynamics of jets in 
quark-gluon plasma. 
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A Energy Loss in a Boosted Fluid 


In this Appendix, we explicitly perform the Lorentz transformation that, as we explained in Sec¬ 
tion 2.2, is needed in order to determine the rate of energy loss dEjdx in the collision center-of- 
mass frame (the frame in which we do the overall computation of the modifications to the energies 
of the partons in the shower) from the rate of energy loss in the local fluid rest frame (the frame 
in which the fluid at the spacetime location of a particular parton is at rest). This transformation is 
particularly important for partons with significant rapidity, as they propagate through fluid that is 
moving with a significant velocity, meaning that the boost from the local fluid rest frame back to 
the collision center-of-mass frame is substantial. 

In the local fluid rest frame, the change in the four-momentum of a parton which propagates 
for an infinitesimal time dtp is given by 

dP^ = Pf {xf, O ^dtF , (A.l) 

where and and Ep are the four momentum and energy of the parton in the local fluid rest frame 
and Ep E^ is the functional form of the rate of energy loss in that frame, in the notation that 
we introduced in Eq. (2.3), and is given by the right-hand side of (2.1) in our hybrid model or by 
the right-hand side of one of the expressions (2.2) in our control models. In writing the expression 
(A.l) we have used our assumption that the exchanges of momentum and energy between the 
parton and the medium do not change the direction of the parton significantly. 

The Lorentz structure of the expression (A.l) simplifies the boost back to the collision center- 
of-mass frame. First of all, it is easy to show that 


dtp dt 
E' 


(A.2) 
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where t and E are the time and energy in the collision center-of-mass frame. Second, after a 
Lorentz transformation, Pj^ and dP^ dP^, with P^ and dP^ the four-momentum and 

the infinitesimal four-momentum loss in the collision center-of-mass frame. Therefore 

^pii pii 

— (A3) 

where the arguments of Pp are still expressed in terms of quantities in the local fluid rest frame. 

We next express the initial parton in the local fluid rest frame as a function of quantities in 
the collision center-of-mass frame via the Lorentz transformation 


E'in = E'in 7F (1 - wv) , (A.4) 

where w = P/E’ is the velocity of the parton in the collision center-of-mass frame and v is the 
local velocity of the fluid in the same frame, which is to say it is the velocity vector for the boost 
between the two relevant frames, is the gamma factor of the local fluid velocity v. 

The relation between xp and x, the distances travelled in the two frames, requires further 
discussion. In the derivation of the rate of energy loss by one of the partons in the shower, xp is 
the distance that that parton has travelled through the fluid. However, as partons propagate through 
the hot plasma created in an heavy ion collision, the temperature and velocity of the fluid at their 
location in space and time changes. We will assume that xp is the accumulated distance of the 
parton summed in such a way that each infinitesimal contribution dxp is evaluated in the local 
fluid rest frame. This means that (if the rate of energy loss depends on xp, as in (2.1) and the 
first expression in (2.2)) the energy lost by a parton traversing some dxp depends on the total 
Xp accumulated by that parton over its previous passage through the flowing plasma. With this 
prescription, which neglects gradient effects, we have 

d:Kp = wdt + 'jp {wp —w)dt^ (A.5) 


where ,l are the transverse and longitudinal components of the parton velocity in the collision 
center-of-mass frame. After some algebra, the increment in the accumulated distance is given by 

^ -h 'jp (v^ — 2vw -h (vw)^) . (A.6) 

Summing over the previous history of the parton, we obtain 

xp{t) = [ dtJ [w^ -h (v2 - 2vw + (vw)2)] , (A.7) 

where to is the creation time of the parton. We have used (A.3), (A.4) and (A.7) in (2.5), (2.6) and 
(2.7). 



B Update on Single-jet and Dijet Observables at ^/s = 2.76 ATeV 

In this Appendix, we update our analysis of the single-jet and dijet observables that we presented 
in Ref. [23] to include the effect of fluid flow on the rate of energy loss, as described in Section 2.2. 
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Figure 13. Jet Raa as a function of for different centralities from our hybrid strong/weak coupling model 
(colored bands) compared to preliminary CMS data from Ref. [103]. The value of the single parameter in 
the model, is fitted to the left-most data point in the top-left panel, namely the jets with 100 GeV < 
Pt < 110 GeV in the most central collisions. All the rest of the features of the colored bands are results 
from our hybrid model. In this Appendix and throughout this paper, single-jet, dijet, photon-jet, Z-jet and 
fragmentation function observables are all fully specified once the single parameter in the model has been 
fixed. 


We also employ the viscous hydrodynamic simulations of Ref. [87], as described in Section 3.1 
and as in all results presented in this paper. As discussed in Section 2.3, we have followed the 
same fitting procedure as we used in Ref. [23] to determine the value of the single parameter which 
controls each the rate of energy loss in our hybrid model and in our two control models. The results 
of these fits are summarized in Table 1. Since the change in all the observables is minor, in this 
Appendix we will not describe each observable in detail. We simply plot the results in Figs. 13, 14, 
15, 16 and 17, and refer the reader to the our previous publication for all discussion. 
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Figure 14. Predictions of our hybrid strongly coupled model for jet i?AA for central Pb-Pb collisions at the 
LHC with ^/s = 2.76 ATeV as a function of px, extended down to px = 15 GeV (left) and for central Au-Au 
collisions at RHIC with ^/s = 200 AGeV (right). In both cases, we only show our results for collisions in 
the 0-10% centrality bin. 
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Figure 15. Dijet imbalance Aj in heavy ion collisions with y/s = 2.76 ATeV from our hybrid model 
(colored bands) compared to CMS data from Ref. [7] (black) and PYTHlA-generated proton-proton collisions 
(violet). Both the hybrid model calculations and the proton-proton reference are smeared according to the 
prescription in Ref. [114]. 


- 41 - 










0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 


1.2 


1 


0.8 
I 0-6 

0.4 


0.2 


0 

100 120 140 160 180 200 220 240 260 280 300 

Pt (GeV) 



30-50% Centrality 
= 2.76 ATeV 





Figure 16. Upper panel: Jet Raa as a function of px for heavy ion collisions with ^/s = 2.76 ATeV in 
two different centrality bins for the hybrid model with its strongly coupled rate of energy loss and the two 
control models from Section 2.1, as compared to preliminary CMS data [103]. Lower panel: Dijet imbalance 
distribution in two different centrality bins for those three energy loss models, as compared to CMS data from 
Ref. [7]. Each of the three models for the rate of energy loss dEjdx includes one free parameter, and in 
each case we have fitted the value of this parameter to obtain agreement between the model and the data for 
100 GeV < px < 110 GeV in the most central (0 — 10%) collisions. 
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Figure 17. Partonic fragmentation function ratios (fragmentation function for jets in Pb-Pb collisions 
over that for jets in p-p collisions) for jets with 100 < < 300 GeV in heavy ion collisions with 

^/s = 2.76 ATeV with four different centralities from the hybrid model and the two control models of the 
rate of energy loss. Jets are reconstructed with the anti-/cT algorithm with R = 0.3. The jet fragments consist 
of final state partons within a cone of radius r = 0.3 around the jet axes determined by the reconstruction 
algorithm. These partons are classified with respect to the longitudinal variable z = pH with pH the 
component of the momentum of the fragment along the jet axis. Note that in the softest region of the 
fragmentation function, say ln(l/ 2 :) > 3.5 — 4, there is an additional contribution to the fragmentation 
function that is not included in the model: as the jet deposits energy and momentum into the medium, some 
part of the medium ends up moving as a wake in the direction of the jet. This will serve to push up the softest 
region of the fragmentation function in all our models, by an amount that should be model independent to 
a good approximation. This effect is not included in our models, but does contribute in the experimental 
data. The three models are compared to the data from Ref. [15]. The separation between the models is 
cause for optimism that the higher statistics measurements expected from LHC heavy ion Run 2 may serve 
to distinguish between models. This optimism must be tempered, however, given that the predictions of the 
models do not differ in a qualitative way and given that we have left out both hadronization and the response 
of the medium to the passage of the jet. 
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Figure 18. Hybrid model prediction for Jet R^a as a function of for different centralities at = 
5.02 ATeV. The single model parameter has been fitted to y/s = 2.76 ATeV data previously and no additional 
parameters have been introduced. 


C Predictions for Single-jet and Dijet Observables at = 5.02 ATeV 

In this Appendix, we provide predictions for the single jet and dijet observables that we presented 
in detail in our previous publication [23] and updated in Appendix B, now for Pb-Pb collisions 
at y/s = 5.02 ATeV. We study 10^ dijet events in p-p collisions at y/s = 5.02 TeV generated 
by Pythia 8.183 [86] without any underlying event. For each of the centrality bins we consider, 
we embed these hard scattering processes into the hydrodynamic simulations along the lines of 
those in Ref. [87] as we described in Section 2.3, but simulations of heavy ion collisions with 
y/s = 5.02 ATeV. We then follow the procedure for determining the energy loss outlined in Section 
3.1. As in our previous computations, we smear our predictions to simulate resolution effects. 
However, and similarly to the procedure we have followed in obtaining our photon-jet and Z- 
jet predictions described in Sections 3.2 and 3.3, we use the smearing functions determined at 
y/s = 2.76 ATeV, since they are yet unknown at ^/s = 5.02 ATeV. 

Fig. 18 shows the hybrid strong/weak coupling model predictions for the suppression factor 
of jets in heavy ion collisions with ^/s = 5.02 ATeV reconstructed with the anti-kp algorithm 
with i? = 0.3. We do not fit even a single parameter here, or anywhere in this Appendix. The 
model is fully constrained by the data on collisions with ^/s = 2.76 ATeV, with its one parameter 
having been fixed as described in Section 2.3 and Appendix B. The width of the displayed bands is a 
combination of our theoretical uncertainties (estimated by varying the temperature Tc below which 
we stop quenching, see Section 2) together with the experimental uncertainties in the data that we 
use to fix the one parameter in each model. Both the centrality dependence and the transverse 
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Figure 19. Hybrid model predictions for the dijet imbalance Aj in Pb-Pb collisions at ^/s=5.02 ATeV 
with different centralities. Both the theoretical calculations and the proton-proton reference are smeared 
according to the prescription in Ref. [114]. 


momentum dependence of in collisions with ^/s = 5.02 ATeV are very similar to what we 
have seen previously at ^/s = 2.76 ATeV, with a slight increase in the suppression of at the 
higher collision energy. Similar conclusions can be drawn from the strong coupling predictions 
for the dijet imbalance in collisions with ^/s = 5.02 ATeV, displayed in Fig. 19. The centrality 
dependence and the momentum dependence of this observable are very similar at the two collision 
energies as well. 

The sensitivity of these predictions to the form we assume for the rate of energy loss is exam¬ 
ined in Fig. 20. As at the lower collision energy, this set of predictions shows little discriminating 
power to our choice among the three models for the rate of energy loss that we have investigated. 
The strongly coupled form (2.1) for dEjdx leads to a slightly large suppression and a slightly big¬ 
ger dijet imbalance than the two control models (2.2), bottom-left panel of Fig. 20, but the effect 
is small compared to current uncertainties. Some separation among models is observed with in¬ 
creasing centrality, top right panel of Fig. 20, although the largest centralities are more and more 
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Figure 20. Model dependence of dijet observables for Pb-Pb collisions with ^/s = 5.02 ATeV. Upper 
panel: Jet R^a as a function of for LHC collisions in two different centrality bins for the three energy loss 
models from Section 2.1. Lower panel: Dijet imbalance distribution in two different centrality bins for the 
three energy loss models. 


sensitive to the energy lost by energetic partons during the hadronic phase of the collision, which 
we are neglecting. 

Finally, in Figure 21 we compare the predictions for partonic fragmentation functions ob¬ 
tained from the three models for the rate of energy loss. These are, again, similar to one an¬ 
other. The modest separation between the model predictions at intermediate values of 2 ; observed 
at ^/s = 2.76 ATeV is also observed at ^/s = 5.02 ATeV. However, the model predictions separate 
most in the region of log z > 3.5 where they cannot be relied upon. In this softest region of the 
fragmentation function, there is an additional contribution that is not included in the model: the 
backreaction of the medium to the jet passing through it will result in additional soft particles in 
the jet cone. 

D Model Dependence of Boson-Jet Correlations 

In this Appendix, we study the sensitivity of the different photon-jet and Z-jet observables that we 
have considered in Sections 3.2 and 3.3 to the microscopic mechanism responsible for energy loss. 
To do so, we repeat our analysis with the strongly coupled form for the rate of energy loss (2.1) 
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Figure 21. Partonic fragmentation functions for jets of 100 < < 300 GeV in heavy ion collisions with 

= 5.02 ATeV for three different models for the rate of energy loss and for four different centralities. Jets 
are reconstructed with the anti-/cT algorithm with R = 0.3. The jet fragments consist of final state partons 
within a cone of angle r = 0.3 around the jet axes determined by the reconstruction algorithm. These partons 
are classified with respect to the longitudinal variable 2 : = pH with pH the momentum of the fragments 
along the jet axis. 


that we use in our hybrid model replaced by one or other of the expressions (2.2) that define our 
two control models, which are inspired by the radiative and collisional energy loss mechanisms. 
As in the strong coupling case, these expressions for dEjdx also contain a single parameter which 
is fitted to the single-jet production rate at one transverse momentum and centrality. The values of 
the parameters obtained from these fits are summarized in Table 1. 

In Fig. 22 we show the results that we have obtained from the hybrid strong/weak coupling 
model and our two control models for the photon-jet observables that we have analyzed in heavy ion 
collisions with ^/s — 2.76 ATeV. The procedure for determining each of the observables is identical 
for all models and is described in Section 3.2. These plots make clear that the uncertainties in the 
present low-statistics data are too large to make it possible to use the photon jet imbalance or the 
spectrum of jets produced in association with isolated photons or the fraction of isolated photons 
with an associated jet to differentiate between microscopic models of the dynanfics of energy loss. 
Within the current uncertainties, all the models agree with the CMS data displayed in these plots. 
The photon-jet imbalance in central collisions, top left panel of Fig. 22, does have modest power 
to discriminate between the strongly coupled model and the control models. This indicates that the 
higher statistics photon-jet data sets anticipated in LHC heavy ion Run 2 could shed light on the 
microscopic dynamics via which jets lose energy. 
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Figure 22. Computations of several photon-jet observables using three different models of the energy loss 
mechanism for heavy ion collisions with ^/s = 2.76 ATeV. The distributions of the transverse momentum 
imbalance of photon-jet pairs for two different centralities are displayed in the upper panels. The middle 
panel shows the ratio of the transverse momentum spectra of jets produced in association with an isolated 
photon in Pb-Pb collisions to that in p-p collisions for two different centralities. The lower panel shows the 
fraction of isolated photons produced in association with a hard jet with > 30 GeV at an azimuthal 
angle more than Ttt/S away from that of the isolated photon. Data are taken from Ref. [111]. 


To better compare our computations with future higher statistics data from LHC heavy ion 
Run 2, we also explore the model predictions of the different energy loss mechanisms for both 
photon-jet observables, displayed in Fig. 23, and Z-jet observables, displayed in Fig 24, in heavy 
ion collisions with ^/s = 5.02 ATeV. As at the lower collision energy, little discriminating power 
is observed. Again as at the lower collision energy, there is some separation among the predictions 
of our hybrid model with its strongly coupled rate of energy loss and the control models in the 
photon-jet and Z-jet momentum imbalance distributions in the most central collisions, displayed in 
the upper-left panels of Figs. 23 and 24, respectively. 
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Figure 23. Predictions for several photon-jet observables computed with three different models of the 
energy loss mechanism in heavy ion collisions with ^/s = 5.02 ATeV. The distributions of the transverse 
momentum imbalance of photon-jet pairs for two different centralities are displayed in the upper panels. 
The middle panel shows the ratio of the transverse momentum spectra of jets produced in association with 
an isolated photon in Pb-Pb collisions to that in p-p collisions for two different centralities. The lower panel 
shows the fraction of isolated photons produced in association with a hard jet with > 30 GeV at an 
azimuthal angle more than Ttt/S away from that of the isolated photon. 
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Figure 24. Predictions for several Z-jet observables computed with three different models of the energy loss 
mechanism in heavy ion collisions with ^/s = 5.02 ATeV. The distributions of the transverse momentum 
imbalance of Z-jet pairs for two different centralities are displayed in the upper panels. The middle panel 
shows the ratio of the transverse momentum spectra of jets produced in association with a Z-boson in Pb-Pb 
collisions to that in p-p collisions for two different centralities. The lower panel shows the fraction of Z- 
bosons produced in association with a hard jet with > 30 GeV at an azimuthal angle more than Ttt/S 
away from that of the isolated photon. 
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